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ABSTRACT 

Many  Navy  applications  of  composite?,,  including  ships'  superstructures, 
submarine  air  flasks  and  missile  rocket  motor  casings,  require  high 
strength  and  reliable  materials.  Composite  strength  reliability  is  dictated  by 
individual  fiber  breaks  at  low  loads  (lower  tail)  and  the  accumulation  of  the 
fiber  failure  sites. 

This  study  examined  the  effects  of  applying  a  preload  to  a  graphite/epoxy 
composite  tow  prior  to  complete  polymerization  of  the  matrix.  The 
objective  was  to  break  the  (inevitable)  weak  fibers  and  minimize  the  effects 
of  the  associated  stress  concentrations,  subsequently  limiting  the  clustering  of 
fiber  failures.  By  eliminating  the  lower  tail,  the  shape  of  the  Weibull 
distribution  is  reduced,  thereby  enhancing  composite  reliability. 
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I.  INTRODUCTION 

In  the  1960s,  Composite  Materials  were  referred  to  as  "Space  Age 
Materials".  Today,  composite  structures  are  common  through  out  our 
society.  The  rapid  growth  of  the  composite  materials  industry  has  been 
phenomenal,  perhaps  only  exceeded  by  the  advancements  in  computer 
technology.  The  concept  of  composite  materials  is  simple,  in  that  two  or 
more  structurally  compatible  materials  are  combined  to  produce  an  end 
product  with  attributes  greater  than  the  sum  of  the  constitutes.  But  as  new 
applications  for  composite  materials  are  developed  so  follows  an  increase  in 
unanswered  questions  with  regards  to  reliability,  especially  associated  with 
vital  national  defense  systems,  public  transportation  systems  and  aerospace 
and  hydrospace  exploration.  As  engineers  it  is  important  that  we  do  not  let 
the  state  of  the  art  exceed  the  state  of  the  science.  We  must  have  answers  to 
those  difficult  questions  concerning  composite  material  applicability, 
maintainability  and  reliability. 

Composite  reliability  is  the  focus  of  attention  at  the  Advanced 
Composites  Laboratory  at  the  Naval  Postgraduate  School,  Monterey 
California.  Here  questions  are  being  investigated  regarding  composite 
reliability  as  a  function  of  service  life,  strength  and  test  methodology.  The 
primary  goal  of  composite  material  research  is  to  be  able  to  have  a  better 
understanding  of  how  composite  materials  behave  over  an  extended  service 
life  and  to  further  investigate  the  mechanics  of  composite  strengthening.  A 


secondary  function  is  the  investigation  into  improved  testing  methods  so  that 
research  techniques  can  be  translated  into  production  applications. 

The  primary  objective  of  this  study  was  to  investigate  composite 
reliability  as  a  function  of  the  strengthening  mechanisms.  This  was 
accomplished  by  studying  the  affects  of  preloading  a  graphite  bundle  prior  to 
complete  curing  of  the  epoxy  matrix.  After  the  matrix  had  cured,  the 
graphite  strand  ultimate  tensile  strength  was  determined  and  the  results  were 
compared  to  a  benchmark  value  based  on  previous  testing. 

In  addition,  individual  fiber  testing  was  conducted  so  that  a 
mathematical  model  may  be  developed  in  order  to  determine  the  optimum 
bundle  preload  level.  Fiber  testing  was  also  conducted  in  order  to 
demonstrate  that  the  Integrated  Fiber  Testing  system  developed  at  the  Naval 
Postgraduate  School  (NPSIT)  could  be  use  as  a  viable  and  accurate  test 
method.  This  was  accomplished  by  conducting  fiber  testing  using  an 
INSTRON  Model  4200  materials  tester  using  current  American  Society  of 
Testing  and  Materials  (ASTM)  procedures  and  comparing  the  results 
obtained  using  the  NPSIT  system.  Graphite  fibers  having  different  strength 
properties  were  tested  in  order  to  show  that  the  testing  systems  used  could 
accurately  ascertain  any  small  differences  in  composite  properties.  The 
focus  of  this  study  was  on  fiber  diameter  and  fiber  tensile  strength. 


II.  BACKGROUND 


Load  bearing  graphite  composite  materials  consist  of  high  strength  /  high 
modulus  graphite  impregnated  with  a  ductile  reinforcing  matrix.  The 
performance  characteristics  of  the  composite  material  are  controlled  by 
three  factors: 

-  Strength  and  modulus  of  the  graphite  fiber 

-  Strength  and  modulus  of  the  epoxy  matrix 

-  Effectiveness  of  the  bond  between  the  matrix  and  fiber,  with 
regards  to  the  mechanics  of  load  transfer  at  the  interface. 

The  overall  tensile  failure  process  of  this  simple  system  is  extremely 
complex,  however  for  the  purpose  of  this  study  this  can  be  simplified  by 
assuming  a  two  dimensional  model  and  that  failure  initiation  occurs  within 
the  fiber  only.  This  failure  model  was  first  presented  by  B.W.  Rosen  (Ref.l 
and  2).  It  consisted  of  parallel  fibers  in  a  homogeneous  matrix,  loaded  in 
tension,  assuming  a  uniform  strain  supported  primarily  by  the  fibers. 

The  most  important  concepts  developed  in  this  classical  study  were  the 
fact  that  the  fibers  were  considered  to  have  statistically  distributed  flaws  or 
imperfections.  Theses  flaws  resulted  in  fiber  failure  at  various  stress  levels, 
therefore  Rosen  considered  each  fiber  as  a  series  of  links  and  that  the 
individual  fiber  failure  process  was  a  function  of  the  weakest  link  principle. 
Finally,  that  the  statistical  strength  distribution  of  theses  flaws  could  be 
approximated  by  a  Weibull  distribution. 
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A.     FIBER  FAILURE  PROCESS 

Initially  when  a  fiber  bundle  ,without  epoxy,   is  loaded  in  tension,  the 

individual  fiber  failures  occur  randomly  at  the  site  of  a  flaw  at  a  load  much 
less  than  failure  load  ( xj  <  Xf  ).  When  this  occurs  the  stress  (ai)  in  the  fiber 

becomes  zero,  and  the  stress  in  the  surviving  fibers  increase  to  a  new  stress 
level  equal  to  some  value   K  *  <j\.   The  existence  of  some  individual  fiber 

failures  may  not  be  catastrophic  in  that  the  increase  in  stress  is  shared  equally 
by  the  remaining  fibers  as  long  as  K  *  01  <  Gf  Then  applied  tension  x  can  be 

increased  to  a  new  value  X2  >  xj  then  this  process  is  repeated  until     K^  *  <j[  > 

of.   Excluding  any  effects  due  to  twist  or  friction,  catostophic  bundle  failure 

occurs  after  each  individual  fiber  breaks  once. 

When  the  fiber  bundle  is  impregnated  with  an  binder  matrix  the 
mechanical  behavior  of  the  composite  strand  is  altered.  The  addition  of  the 
matrix  creates  a  unique  load  transferring  system.  This  provides  for 
transverse  load  sharing  between  each  individual  fiber  and  longitudinal  load 
sharing  to  another  segment  of  the  fiber.  The  matrix  also  provides  the 
mechanism  for  localizing  the  effects  of  microcracks  within  the  composite 
material.  These  two  functions  combined  provide  for  the  distinct  mechanical 
performance  of  fiber  reinforce  composite  materials. 

The  composite  failure  process  based  on  Rosen's  failure  model,  with 
matrix  binder,  predicts  different  results  from  dry  bundles.  When  the  initial 
load  Xi  is  applied  initial  fiber  failures  still  occur  at  the  inevitable  flaw  sites, 

however  the  effects  are  confined  to  a  localized  region.  Stress  a  decreases  to 
zero,  and  the  shear  x  in  the  matrix  reaches  a  maximum  value.  The  axial  load 
is  transmitted  by  this  shear  to  only  the  neighboring  fiber  segments,  (see 
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Figure  1).  The  neighboring  fibers  now  show  a  load  concentration  due  to  the 
fiber  break.  At  the  broken  fiber  tip  there  exists  a  length  of  fiber  which  is 
ineffective  in  carrying  the  applied  load.  This  is  known  as  the  ineffective 
length  5.  The  increase  in  shared  stress  on  the  immediate  neighbors  is 
substantially  greater  then  that  shared  by  the  fibers  that  are  a  greater  distance 
from  the  failure  site. 
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Figure  1 .  Fiber  Bundle  Model 
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The  determination  of  the  ineffective  length  8  is  difficult,  and  requires  an 
understanding  of  the  shear  stress  distribution  along  the  interface.  The 
ineffective  length  is  estimated  by  Rosen  [Ref.  1]  to  be: 

5-1/2   {(Vf-l/2-l)Ef/Gm}l/2  cosh-l{[l  +  (l+<p)2]/2[l-q>]}df 

Vf  -  Volume  fraction  of  fiber 

Ef  -  Modulus  of  fiber 

Gf  -  Shear  modulus  of  matrix 

cp  -  Fiber  efficiency  or  the  fraction  of  the  undistributed  stress 

value  below  which  the  fiber  is  considered  to  be  ineffective 
df  -  Fiber  diameter 

B .     COMPOSITE  FAILURE 

In  general  there  are  three  basic  modes  of  composite  failure.  The  first  is 
debonding  which  is  failure  at  the  fiber/matrix  interface  which  is  caused  by  a 
high  interface  shear  stress  (Figure  2a).  This  failure  may  start  with  a  single 
broken  fiber  and  propagates  within  the  interface,  along  the  fiber  length, 
which  in  turn  drastically  increases  the  ineffective  length  5.  Thus  increasing 
the  likelyhood  of  the  crack  propagating  to  another  flaw  site.   The  study  of 

this  type  of  failure  identifies  the  need  for  improving  the  interface  properties 
such  as  adhesion  and  toughness  properties  of  the  epoxy  (Note:  Xj  «  Xf).  The 

second  type  of  composite  failure  is  crack  propagation  (Figure  2b).  This 
occurs  when  an  initial  crack  is  developed  at  the  flaw  site  and  propagates 
transversely  across  the  composite.  Crack  propagation  is  controlled  by  the 
fracture  toughness  of  the  matrix  and  fiber  and  thus  minimizing  the  effects  on 
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local  load  sharing.  This  failure  mode  may  be  associated  with  some  slight 
debonding.  (Note:  X[ «  Xf). 


f"N   s~^ 

p^ 

Figure  2a. 


Figure  2b. 


Figure  2c. 


Figure  2.  Composite  Failure  Model 

Finally,  if  debonding  and  crack  propagation  are  controlled  then 

composite  failure  occurs  when  as  load  is  increased,  the  number  of  statistically 

distributed  fiber  failures  accumulates  in  the  vicinity  of  one  cross  section  such 

that  when  the  load  shared  in  the  surviving  fibers  exceeds  the  ultimate  load  of 
the  system  (  X[  >  Xf ).    This  mode  of  failure  is  usually  a  combination  of 

debonding  and  crack  propagation,  occurring  at  a  much  higher  load  level, 
(Figure  2c). 
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C.      LOAD  SHARING 

The  most  pronounced  effects  of  adding  a  matrix  binder  to  the  fiber 
bundle  is  that  of  load  sharing.  This  is  best  discussed  using  a  microbundle 
model  developed  by  S.L.  Phoenix.  [Ref.  3].  This  model,  (Figure  3)  is  made 
up  of  n  fibers,  each  consisting  of  m  links  having  a  length  of  5,  the  ineffective 
length.    In  each  microbundle  (  =    n  *  5)  there  exists  a  set  of  statistically 

distributed  flaws.  When  the  load  x  is  applied  the  fibers  will  sequentially  fail, 
in  accordance  with  a  Weibull  distribution  function. 

F(x)  =  1  -  exp  [-(x/(3x)«] 

Where  a  is  the  shape  parameter  and  (3  is  the  scale  parameter. 

Because  the  individual  fiber  consists  of  m  links  ,  the  weakest  link 
theorem  applies.  Therefore,  the  value  of  (3X  is  a  function  of  gage  length  and 

must  be  corrected  for  GL. 

P2x    =    Plx    (GLi/GL2)l/<* 

When  an  initial  fiber  failure  x\  occurs  in  a  weakest  link  the  load  on  each 
immediate  neighbor  increases  to  a  value  X=Ki*x,  where  K\  is  the  load 
concentration  factor  and  Kj>1,  (Figure  3).  If  another  weak  link  exists  within 

the  region  of  concentrated  load,  then  a  second  fiber  failure  will  prematurely 

occur.     This  in  turn  creates  an  additional  load  sharing  burden  on  the 
neighboring  fibers    X  =  K2*x,  where  K2>Ki>l,  (Figure  3).    Therefore, 

randomly  created  load  concentrations  may  prematurely  break    a  larger 
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number  of  fibers  at  a  load  x\  than  would  otherwise  normally   occur, 

subsequently  weakening  the  composite.   J.M.  Hedgepeth  [Ref.  4]  provided 
quantitative  values  to  the  load  concentration  factors,  where  K\  =  1.33,  K2  = 

1.60,  K3  =  1.83  and  K4  =  2.03.  Meaning,  for  example,  that  the  load  in  the 
fibers  adjacent  to  three  consecutive  fiber  breaks  equals  1.83  *  applied  load  . 
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Figure  3.  Load  Concentrations  Model 

If  the  random  fiber  breaks  occuring  at  X[  were  dispersed  throughout  the 

composite  the  effects  of  load  concentrations  K  would  be  minimal,  however 
this  can  not  be  controlled.  Therefore  an  attempt  must  be  made  to  control  the 
effects  of  the  spatial  clustering  of  fiber  breaks  and  load  concentrations  in  a 
local  region. 
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A.S.  Tetelman  [Ref.  5]  and  C.  Zweben  [Ref.  6]  provide  excellent 
summaries  of  the  composite  failure  process  and  the  effects  of  load 
concentrations.  Tetelman  also  discussed  a  few  concepts  for  minimizing  load 
concentrations  and  controlling  crack  propagation.  These  include  the  use  of 
discontinuous  fibers,  using  fibers  or  matrix  consisting  of  various  modulus  of 
elasticity  and  the  use  of  prestressing  at  elevated  temperatures  In  fact 
prestressing  composite  bundles  at  elevated  temperatures  was  proven  to  be  a 
viable  solution  during  experimental  work  conducted  by  E.M.  Wu.  [Ref.  7]. 
GJ.  Mills  [Ref.  8]  has  also  studied  the  effects  of  prestressing  Boron/Epoxy 
prepregs  by  rolling. 

D.     COMPOSITE  RELIABILITY 

The  reliability  of  a  composite  structure  depends  on  a  definite 
characterization  of  the  weak  strength  distribution  (the  weak  lower  tail). 
Definitive  evaluation  of  the  lower  tail  requires  a  large  sample  size  and 
consistent  material  production  system. 

Such  a  data  base  for  composite  material  performance  is  usually 
developed  from  either  limited  fiber/bundle  testing  or  complete  system  testing 
of  a  small  sample  size.  A  small  data  base,  when  combined  with  the  existence 
of  the  statistically  distributed  flaws  or  imperfections  in  the  graphite  fiber 
due  to  manufacturing,  processing  or  handling  practices,  makes  the  behavior 
of  a  given  composite  structure  difficult  to  predict.  When  a  multimillion 
dollar  rocket  motor  casing  or  aircraft  component  consists  of  an  enormous 
number  of  fiber  links,  performance  reliability  becomes  dependent  upon  a 
statistical  study  of  the  strength  characteristics.    Because  of  the  existence  of 
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the  fiber  flaws,  the  strength  values  of  the  Weibull  distribution  are  well 
dispersed  within  the  lower  tail.  This  imposes  great  restrictions  on  the 
allowable  design  limits  for  the  composite  structure.  Ideally,  there  would 
exist  nearly  no  dispersion  in  the  test  data  (coefficient  of  variance  approaches 
zero)  and  then  the  data  would  be  completely  reliable  for  a  given  strength 
level. 

This  study  will  focus  on  improving  composite  strength  reliability  by  the 
use  of  preloading  the  fiber  bundle  prior  to  the  curing  of  the  epoxy  resin  in 
order  to  remove  the  weak  fiber  sites  (lower  tail  of  the  Weibull  distribution) 
and  decrease  the  dispersion  of  the  composite  strength. 
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HI.  EXPERIMENTATION 

The  experimental  work  was  concentrated  in  two  areas.  First, 
determining  the  diameter  and  strength  of  the  individual  fiber  and  second, 
determining  the  strength  characteristics  of  a  graphite  tow  subsequent  to 
preloading. 

The  material  tested  was  from  a  Hercules  Magnamite  high  strength 
graphite,  type  AS-4  spools  008  and  019.  The  graphite  bundle  consisted  of 
3000  fibers  with  a  nominal  diameter  of  7  micrometers  and  a  denier  of 
.005746  grams/inch. 

A.     FIBER  TESTING 

Individual  fiber  testing  was  a  vital  part  of  this  study  which  provided  the 
essential  data  required  to  determine  the  appropriate  preload.  First  the  fiber 
diameter  was  determined  then  the  fiber  failure  strength.  This  information 
was  used  to  develop  a  Weibull  probability  plot  to  statistically  determine  the 
number  of  fiber  failures  as  a  function  of  applied  load. 

To  study  the  feasibility  of  using  an  alternative  testing  method,  two 
testing  methods  were  used  during  this  study.  The  first  method  was  the  use  of 
an  Instron  Material  Tester  (INSTRON),  (Figure  4),  with  associated  ASTM 
procedures  to  determine  fiber  failure  load.  Second,  using  a  Integrated 
Testing  system  under  development  at  the  Naval  Postgraduate  School, 
Monterey,  California.  (NPSIT),  (Figure  5). 
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Figure  4.  INSTRON  Testing  System 


Figure  5.    NPSIT  Testing  System 
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The  primary  advantage  to  the  NPSIT  system  is  that  the  fiber  diameter  as 
well  as  the  failure  load  can  be  readily  determined  on  the  same  test  stand.  This 
minimizes  the  handling  of  the  test  samples  and  bias  the  statistical  parameters. 

The  correlation  between  fiber  test  length  and  failure  load  has  been  well 
established  in  past  studies  and  is  not  addressed  here.  A  gage  length  of  5  centi- 
meters was  used  for  fiber  testing.  Strain  rate  or  crosshead  speed  also  has  an 
effect  on  test  results,  however  this  effect  was  not  studied.  A  constant  speed  of 
10  percent  of  gage  length  per  minute  (in  accordance  with  ASTM  procedure) 
was  employed. 

Regardless  of  the  testing  method  used,  sample  preparation  (Appendix 
A.)  and  testing  parameters  were  identical. 

1.       Fiber  Diameter 

Current  methods  for  determining  fiber  diameter,  as  outlined  in 
ASTM  procedures,  call  for  the  use  of  the  Scanning  Electron  Microscope 
(SEM).  This  is  a  time  consuming  and  some  what  arduous  task  when  testing 
large  numbers  of  fibers. 

Recent  work  at  the  Naval  Postgraduate  School,  Monterey,  has  provided  two 
additional  methods  which  are  more  conducive  to  this  study.  The  principles 
and  procedures  as  developed  by  T.  A.  Bennett  [Ref.  9]  are  based  on  the 
physical  optics,  that  if  a  illuminated  beam  of  light  is  directed  at  the  fiber  that 
a  diffracted  light  pattern  is  produced  which  is  perpendicular  to  the  fiber 
( Figure  6). 
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Figure  6.  Diffraction  Pattern 

The  first  method  Bennett  looked  at  was  the  use  of  a 
photoconductive  cell  to  measure  the  distance  between  nodes  (null  intensity  in 
the  diffracted  light  pattern).  The  second  method  consists  of  the  use  of  a 
microneye  interfaced  with  an  Apple  Plus  II  microcomputer.  Bennett  [Ref.  9] 
and  M.  Storch  [Ref.  10]  then  developed  the  data  acquisition  techniques 
necessary  to  accurately  determine  fiber  diameter.  The  required  test 
configuration  is  essentially  the  same  for  both  methods,  however,  since  only 
the  photocells  was  used  in  this  study  ,  only  those  associated  procedures  will  be 
further  discussed. 

The  test  equipment  (  Figure  5)  consists  of ;  a  low  power  (.052  mW) 
Helium-Neon  laser,  which  provides  the  required  light  source  of  known  wave 
length  ( lambda  =  632.8EE-09  meters).  A  focusing  lens,  a  spatial  filter  and  a 
collimating  lens  are  enclosed  as  one  unit  and  mounted  on  the  end  of  the  laser. 
This  is  necessary  in  order  to  produce  the  parallel  light  pattern.  The  fiber  is 
mounted  on  the  test  stand    which  by  use  of    micrometers  allows  for 
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longitudinal  and  transverse  adjustments,  rotation  of  the  sample  and 
adaptability  to  the  load  testing  device.    At  the  end  of  the  tracks  are  two 
photocells  mounted  to  adjustable  pedestals  which  allow  for  transverse 
adjustments.  The  photocells  are  connected  to  two  Fluke  8840A  multimeters. 
2.      INSTRON  Load  Test 

Individual  fiber  load  testing  was  performed  using  the  Instron  Load 
Tester  Model  4206  (INSTRON)  inconjunction  with  compatible  data 
acquisition  systems.  The  data  acquisition  system  to  be  used  during  this  study 
is  the  Hewlett-Packard  3497A  Data  Acquisition/Control  Unit  and  the 
Hewlett-Packard  HP-85  for  computer  programming.  Appendix  B  provides 
a  detailed  listing  of  the  Interactive  Data  Acquisition  Software  (IDAS).  IDAS 
was  design  to  minimize  operator  effort  and  ensure  accurate  and 
reproducible  test  results.  All  IDAS  procedures  were  written  in  accordance 
with  American  Society  of  Testing  and  Materials  (ASTM)  Standard 
Specifications. 

After  the  required  program  is  loaded  enter  RUN  and  then  follow 
the  outlined  steps  in  Appendix  B.  The  programs  use  the  function  keys  to 
provide  for  flexible  data  acquisition.  Options  are  normally  provided  for  a 
CRT/printer  output  or  with  the  graphics  plotter. 

Fiber  diameter  was  measured  prior  to  this  test  with  the  use  of  the 
NPSIT  test  system. 

The  IDAS  software  includes  the  following  programs. 

LDCALB  -  Load  Cell  Calibration:  Characterization  of  the 
strength  of  a  fiber  required  absolute  calibration  of  the  load  transducer.  This 
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was  accomplished  by  calibrating  the  digital  output  of  the  load  transducer  by 
standard  weights  at  six  levels.  This  provided  absolute  calibration  all 
analogue  and  the  analogue  to  digit  interface. 

MCTST  -  Characterization  of  the  stiffness  of  a  fiber  requied  the 
calibration  of  the  system  compliance.  The  total  compliance  measured  (Ct)  is 
the  sum  of  the  fiber  compliance  (Cf)  plus  the  system  compliance  (Cs),  all 

three  are  a  function  of  applied  load  level  P. 

Ct(p)  =  CKp)  +  Cs(p) 

The  system  compliance  Cs  is  comprised  of  the  compliance  of  the 
testing  machine,  of  the  cardboard  tab  and  the  adhesive  which  bonds  the  fiber 

to  the  tab.  If  the  system  compliance  is  assumed  to  be  constant,  then  a  single 
calibration  constant  (  effective  gauge  length  of  zero)  is  determined  as 
recommended  by  ASTM  procedure  D-3379-75(82),  (Ref.  11).  as  discussed 
in  ASTM  procedure  3379-75(82),  the  effective  gauge  length  can  be 
extrapolated  from  measurement  of  the  total  compliance  using  three  different 
gauge  lengths.  The  short  coming  of  this  procedure  is  that  it  assumes  that 
system  compliance  is  independent  of  load  level  P.  The  latter  was  observed  to 
be  a  gross  assumption.  Compliance  calibration  was  accomplished  using  an 
actual  zero  gage  length  sample.   In  this  case,  the  fiber  compliance  Cf  term 

vanishes  and  the  total  compliance  measured  is  the  system  compliance. 
Furthermore,  the  compliance,  over  the  entire  load  range  of  the  strength  of 
the  fiber  measured  was  recorded  and  the  data  fitted  by  a  second  degree 
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polynomial.  This  second  order  polynomial  was  then  used  to  calibrate  the 
tensile  strength  test  data,  thereby  compensating  for  the  non-linear  depencey 
on  the  load  level. 

QFIT  -  Quadratic  Curve  Fit:  This  program  reads  the  saved 
Machine  Compliance  data  and  provides  the  three  coefficients  of  a  quadratic 
curve.  These  values  E,  F,  G  were  then  recorded  for  further  use  in  the  fiber 
test  program  to  mathematically  subtract  from  the  total  displacement  the 
displacement  associated  with  the  system  compliance. 

FTST  -  Fiber  Test:  Is  the  primary  program  in  this  series. 
INPUTS,  allows  for  the  input  of  all  key  parameters  necessary  to  determine 
the  actual  fiber  ultimate  strength.  The  test  is  run  for  30  seconds  at  a  speed  of 
.lx  gauge  length  (5  millimeters/minute  on  a  50  millimeter  sample).  The  A, 
B,  C,  E  and  F  inputs  as  discussed  above  are  required  to  perform  the  necessary 
mathematical  calculations.  AQUIR,  provides  for  the  actual  proof  test  and 
data  acquisition.  Once  the  sample  is  ready  the  program  will  adjust  for  the 
initial  tension  then  start  the  timer,  providing  for  300  data  readings  during  the 
30  second,  run.  The  results  can  be  displayed  on  the  CRT/printer  using  GRAF 
or  on  the  HP  Plotter  using  PLOT.  SAV-DAT  and  RD-DAT  provide  for  the 
storage  and  reading  of  test  data.  CAL,  mathematically  subtracts  system 
compliance  from  the  total  displacement  and  provides  the  unbaised  results. 
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LDTST  -  Load  Test:  Provides  for  a  quick  check  of  the  load  cell 
calibration  using  the  original  A,  B  and  excitation  voltage  from  the  LDCALB 
results. 

LDSTOR  -  Load/Store:  is  a  data  management  tool  for  transfering 
files. 

3.      NPSIT  Load  Testing   System 

The  Naval  Postgraduate  School  Integrated  Load  testing  system  was 
developed  as  an  alternative  testing  method  which  provides  for  the 
determination  of  fiber  diameter  as  well  as  failure  load  with  one  test  system. 

The  data  acquisition  and  software  is  essentially  the  same  as  that  used 
with  the  INSTRON  test  system.  The  only  exception  being  that  a  short 
program  was  added  to  FTST  so  that  fiber  diameter  and  failure  load  are 
recorded  on  the  same  file  for  easier  access. 

By  following  the  procedures  in  Appendix  A  and  Appendix  B,  the 
NPSIT  system  can  be  quickly  initialized  and  ready  for  fiber  testing.  After 
the  FTST  program  is  loaded,  mount  the  fiber  sample  in  the  clamps.  First 
measure  the  fiber  diameter,  as  discussed  and  then  enter  AQUIR.  The  distinct 
difference  between  this  test  method  and  the  INSTRON  test  method  is  that  the 
NPSIT  drive  motor  is  manually  controlled  by  the  operator  upon  prompt 
from  the  HP-85  program,  where  as  the  INSTRON  tester  is  software 
controlled. 

All  fiber  data  is  recorded,  displayed  and  stored  in  the  same  manner 
as  with  the  INSTRON  test. 
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B.     TOW  TESTING 

The  objective  of  graphite  tow  testing  was  to  determine  the  the  effects  on 
the  strength  characteristics  as  a  result  of  applying  a  preload.  Experimentally, 
AS-4-019  graphite  bundles,  from  the  same  spool  used  in  fiber  testing,  were 
tested. 

1.  Initial  Set-up 

The  samples  were  prepared  by  paying  the  graphite  out  on  a  table 
and  applying  a  2  kilogram  load.  Then  copper  tabs  (1  inch  x  1.5  inch)  were 
fasten  to  the  bundle  at  a  gage  lengths  of  10  inches  using  an  epoxy  mix  as  the 
adhesive.  The  2  kilogram  load  removes  the  slack  and  aligns  the  individual 
fibers  within  the  bundle.  The  copper  tabs  provide  for  an  area  for  the  bundle 
to  be  gripped  in  the  Instron  tester. 

The  Instron  Universal  Testing  Instrument  Model  4206  was  used  to 
provide  the  tensile  load  to  the  test  samples.  This  system  was  used  in 
conjunction  with  the  4200  Series  Expanded  Control  Console  with  data 
acquisition  being  provided  by  the  use  of  an  IBM  PC- AT.  The  Instron/IBM 
data  acquisition  on  system  provides  for  the  real  time  graphical  display  of  load 
versus  displacement  and  the  creation  of  a  data  base  for  additional 
computations  and  graphical  outputs.  The  initial  operator  and  report  files  are 
established  prior  to  the  commencement  of  testing  to  generate  the  desired 
reports  and  graphs. 

2.  Preloading  Bundles 

Prior  to  applying  a  preload,  the  Graphite  bundle  was  impregnated 
with  the  epoxy  matrix.  The  matrix  consisted  of  a  mixture  of  Dow  DER-332 
resin  (55  percent  by  weight)   and  Texaco  Jeffamine   T403  (45  percent  by 
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The  final  resin  content  of  the  strand  was  determined  by  weighing  a 
few  random  samples  and  comparing  to  the  weight  of  the  dry  bundle. 

%R  =  WT-WP 
WT 

Where    Wj  is  the  total  weight  of  the  strand  and  Wg  is  the  weight  of  the 

bundle  (denier  x  gage  length).  A  resin  content  of  55  to  60  percent  is  desired. 

3.      Load  Testing 

The  failure  load  of  the  final  product  (the  matrix  impregenated 

composite  strand)  was  determined  using  the  Instron/IBM  system  described 

earlier.  After  the  strand  is  mounted  in  the  test  grips  the  desired  test  files  are 

inputted  and  the  machine  crosshead  speed  is  set  a  10  centimeters  per  minute. 

The  test  is  allowed  to  continued  until  complete  tensile  failure  of  the 

composite  strand  occurs.  The  final  output  provides  for  the  maximum  load  in 

kilograms  and  displacement  at  failure  in  millimeters.    A  graphical  report 
may  also  be  generated. 
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VI.  RESULTS 

A.     FIBER  TESTING 
1.      Fiber  Diameter 

A  review  of  the  statistical  summary  of  the  fiber  diameter  results 
(Table  I),  provides  a  clear  indication  that  the  test  method  used  for 
determining  fiber  diameter  produced  consistent  and  accurate  results.  The 
normal  probability  plot,  Figure  8,  shows  how  the  AS-4-008  series  and  the 
AS -4-0 19  series  graphite  diameters  compared  to  the  merged  results. 

A  total  of  79  samples  were  measured  from  the  AS-4-008  graphite 
spool,  with  a  mean  diameter  of  7.233  microns  versus  82  fiber  samples  from 
the  AS-4-019  spool  with  a  mean  diameter  of  7.245  microns.  The  large 
variability  of  fiber  diameter  (as  measured  by  the  magnitude  of  the  standard 
deviation)  verifies  that  the  fiber  diameter  is  not  consistent  over  the  entire 
length. 

Exploratory  effort  was  expended  in  the  adjustment  of  the  test  stand, 
optimizing  the  distance  between  the  fiber  and  the  projected  plane  (the 
photocell)..  This  was  required  to  compensate  for  the  narrow  width  of  the 
photocell  window.  It  was  necessary  to  adjust  the  test  stand  so  the  width  of  the 
node  was  nearly  as  wide  as  the  photocell.  If  the  node  is  too  wide,  there  exists 
several  locations  of  maximum  resistance  readings.  If  the  the  node  is  too 
narrow,  too  much  light  saturates  the  photocell,  thus  making  it  impossible  to 
locate  the  center  of  the  node. 
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TABLE  I:  STATISTICAL  SUMMARY  OF  SINGLE  FIBER  TESTING 
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Figure  8.     Normal  Distribution  Plot  of  Fiber  Diameter 
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2.      Fiber  Failure  Load 

Interpretation  of  fiber  failure  data  is  focused  on  two  areas 
--    Comparison  of  test  results  from  INSTRON  and  NPSIT  testing  to 

identify  any  differences. 
—    Comparison  of  the  statistical  strength  for  samples  from  two  different 
spools. 

In  Table  II  statistical  parameters  from  Table  I  are  employed  for 
comparison  of  the  INSTRON  and  NPSIT  test  methods. (see  Figures  9  and  10). 

TABLE  E.  FIBER  LOAD  RESULTS 
Mean(gmf)  Std  Dev 

008  019  008  019 

INSTRON         14.26  13.91  4.20  4.82 

NPSIT  14.56  14.06  2.78  3.20 

Difference        +2%  +1%  -15%  -15% 

Comparing  the  first  two  columns  of  TABLE  II,  it  is  apparent  that  the  mean 
strengths  of  the  samples  measured  with  the  NPSIT  are,  slightly  but 
consistently,  higher  for  both  spools  tested.  Comparing  the  last  two  columns, 
the  variabilities  of  the  measured  samples  in  the  NPSIT  are,  significantly  and 
consistently,  lower  for  both  spools.  These  observation  substantiate  the 
hypotheses  that  the  NPSIT,  which  integrates  the  fiber  diameter  measurement 
with  the  strength  measurement,  enables  the  elimination  of  several  handling 
and  storage  steps,  thereby  minimizing  the  associated  accidental  damage. 
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Figure  10.     Weibull  Distribution  Plot  of  Fiber  Load 
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Minimization  of  this  damage  results  in  a  shift  of  the  lower  tail,  of  the  Weibull 

plot,  to  the  left  (evident  in  both  Figure  9  and  Figure  10  for  samples  labeled 
*)  and  a  significant  reduction  of  variability  and  the  standard  deviation. 

Minimization  of  the  handling  damage  cannot  increase  the  instrinsic  strength, 
therefore  has  no  effect  on  the  upper  tail  (again  evident  in  Figures  9  and  10). 
The  shift  of  the  lower  tail  changes  the  geometric  centroid,  of  the  probability 
distribution  function  (pdf),  therefore  causing  a  slight  increase  of  the  mean. 
The  consistency  between  the  quantitative  measurement  with  the  qualitative 
expectations  support  the  conclusion  that  the  integrated  NPSIT  test  is  a  worth 
while  improvement. 

For  comparison  of  the  statistical  strength  of  the  two  spools,  the 
strength  is  modeled  by  the  two  parameter  Weibull  distribution,  Table  III. 
The  parameters,  as  based  on  Maximum  Likelihood  Estimation,  of  the 
respective  data  sets  are  summarized: 

TABLE  HI.    FIBER  LOAD  WEIBULL  PARAMETERS 


Alpha 

a 

Beta  p 

NPSIT 

INSTRON 

MERGE 

NPSIT 

INSTRON 

MERGE 

008 

5.75 

3.62 

4.28 

15.72 

15.80 

15.79 

019 

4.86 

3.32 

3.94 

15.33 

15.52 

15.45 

Diff 

-18% 

-9% 

-8.5% 

-2.5% 

-1.8% 

-2.2% 

Again,  note  the  significant  difference  between  the  a  parameter 
(which  is  inversely  associated  with  variability)  for  the  two  spools.  This  is 
particularly  evident  when  handling  has  been  reduced    On  the  other  hand, 
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consistent  with  the  above  discussion,  the  difference  between  the  (3  parameter 
(associated  with  mean  strength)  is  slight.  Physically,  composite  samples 
fabricated  with  the  019  spool  which  has  a  smaller  a  (or  larger  variability) 
will  have  a  larger  number  of  fiber  breaks  than  those  fabricated  with  the  008 
spool.  In  accordance  with  the  local  load  sharing  model,  the  larger  number  of 
fiber  breaks  will  increase  the  stochastic  process  of  clustering  of  failure  sites 
hence,  leading  to  a  higher  composite  scatter.  Therefore,  preloading  spool 
019  and  observed  accompanied  change  of  composite  variability  will  offer 
both  a  verification  of  the  load  sharing  model  and  explore  the  practicality  of 
improving  composite  reliability. 

Complete  tabulated  results  for  all  fiber  diameter  and  tensile 
strength  measurements  are  found  in  Appendix  D. 

B.     TOW  BUNDLE  PRELOAD 

The  objective  of  preloading  the  tow  bundle,  prior  to  matrix  binder 
impregnation  was  to  decrease  the  variability  of  the  composite  strength.  It 
was  established  in  [Ref.l]  and  [Ref.  2]  that  composite  failure  occurred 
sequentially  when  weak  fiber  failure  began  at  low  load  levels.  As  applied 
load  is  increased,  the  number  of  broken  sites  increased  and  formed  spatially 
clustered  failure  sites  leading  to  stress  concentrations,  the  most  severe  of 
which  causes  catastrophic  failure.  This  failure  process  can  be  observed  by 
comparing  the  load  deformation  curve  of  a  matrix-free  tow  bundle  (curve  B 
in  Figure  11)  and  a  matrix  impregnated  bundle,  the  composite  (curve  C). 

For  the  tow  bundle  curve  B,  the  initial  slope  of  the  load  deformation 
curve  is  the  sum  of  all  the  filaments  in  the  bundle.    As  the  deformation 
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Figure   11.     Composite  Deformation  Curve 
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displacement  is  increased,  the  weakest  fibers  break  and  the  load  carrying 
capability  of  the  entire  bundle  is  lost,  thus  resulting  in  a  decrease  in  slope 
such  that  the  secant  modulus  is  equal  to  the  sum  of  modulus  of  the  surviving 
fibers.  This  continues  until  only  a  few  very  strong  fiber  survive  and  the 
secant  of  the  load  deformation  curve  of  the  bundle  approaches  horizontal. 

For  the  composite  curve  C,  the  initial  slope  of  the  load  deformation 
curve  is  the  same  as  that  of  the  bundle  curve  B,  since  both  have  the  same 
number  of  filaments.  Within  the  composite,  weak  fibers  start  to  fail  at  0.3 
percent  strain  as  in  the  bundle.  However,  in  the  presence  of  the  matrix 
binder,  the  load  carried  by  the  broken  fiber  is  transferred  to  the  neighboring 
fibers  via  the  matrix  and,  instead  of  the  entire  fiber,  only  a  portion  equal  to 
the  ineffective  length  is  lost.  The  ineffective  length  is  approximately  equal  to 
10  fiber  diameters  or  in  the  case  of  graphite,  in  the  order  of  lOOum.  Such  a 
minute  loss  can  not  be  detected  from  the  resolution  of  the  load  deformation 
curve.  As  result  the  load-deformation  curve  remains  apparently  linear  up  to 
the  point  of  catastrophic  rupture  even  though  numerous  filament  failures 
sites  accumulated  internally. 

Therefore,  it  is  desirable  to  develop  a  preload  procedure  to  break  the 
inevitable  weak  fibers,  without  creating  stress  concentrations  which  lead  to 
premature  undue  breakage  of  the  neighboring  fibers.  This  is  possible  if  the 
preload  can  be  performed  before  the  matrix  can  effectively  transfer  the  load. 
Subsequently  when  the  matrix  becomes  effective,  the  load  around  the  broken 
site  can  be  transferred  to  the  neighbor  without  causing  spatial  clustering  of 
broken  fiber  sites  and  subsequently  reducing  the  variability  of  the  ultimate 
strength. 
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The  most  obvious  implementation  is  to  perform  the  preload  on  a  dry 
bundle  prior  to  the  introduction  of  the  matrix  binder.    Exploratton  proved 
that  this  approach  was  not  practical.    In  the  absence  of  matrix  bmder,  a 
broken  filament  unleashed  the  stored  elastic  energy  dynamically  leadmg  to 
fraying  and  lost  of  alignment  of  the  filaments,  (see  Figure  12).  In  tact  this 
phenomenon  occurred  when  preload  was  applied  to  a  bundle  wetted  ,n  a 
matrix  but  prior  to  polymerization.  In  such  a  case  the  surface  tenston  of  the 
uncured  polymenc  matnx  is  insufficient  to  contam  the  dynarmc  energy 
release  of  the  breaking  filaments.    Therefore,  the  impregnated  bundle  was 
allowed  to  air  cure  prior  to  preloading  to  ensure  sufficient  surface  tensron  to 
prevent  fraying. 


Figure  12.  Fraying  of  Preloaded  Strands 
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Based  on  the  results  from  Appendix  C,  the  preload  level  at  a 
deformation  of  3.0  millimeters  would  break  40  percent  of  the  fibers.  It  can 
be  expected  that  this  is  the  effective  strain  level  desired  to  verify  the  preload 
mechanism.  Because  of  the  fraying  problem  associated  with  this  method,  the 
preload  level  was  reduced  to  a  level  below  which  fraying  occurred.  After  a 
series  of  samples  were  tested  at  various  strains,  it  was  determined  that  2.25 
millimeters  provided  the  best  results  for  a  range  of  room  temperature  curing 
times  of  19  to  22  hours.  At  2.25  millimeters,  it  can  be  predicted  that  13 
percent  of  the  fibers  were  deliberately  broken. 

1.  Procedure 

The  finalized  preload  procedure  consisted  of  : 

-  Homogeneous  mixing  of  epoxy  and  curing  agent 

-  Impregnate  graphite  tow  with  epoxy  in  a  flat  pan 

-  Partially  cure  epoxy  matrix  at  room  temperature  for  19  hours 

-  Preload  samples,  under  displacement  control  to  2.25mm  ( .9% 
strain) 

-  Final  cure  samples  at  60°C  for  16  hours 

-  Tensile  test  samples  at  crosshead  displacement  of  10  mm/min. 

Ultimate  tensile  test  results  for  the  preloaded  and  non-preloaded 
samples  are  listed  in  Appendix  D.  Table  XV. 

2.  Interpretation 

The  tensile  strengths  associated  with  the  preloaded  sample  are 
compared  with  the  bench  mark  non-preloaded  samples.  The  non-preloaded 
samples  were  prepared  essentially  by  the  same  procedures  listed  above  with 
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the  exception  of  air  curing  and  preloading,  using  an  automated  filament 
winding  set-up  with  no  manual  handling.  The  expected  quality  of  the  non- 
preloaded  samples  is  comparable  to  well  controlled  production  quality.  The 
preloaded  samples  were  prepared  individually  with  all  manual  handling.  The 
quality  was  expected  to  be  lower  than  a  fully  developed  production  process. 
One  can  expect  that  any  improvement  observed  with  the  preloaded  sample 
may  be  realized  in  a  fully  developed  process. 

The  tensile  strengths  of  both  the  non-preloaded  and  the  preloaded 
samples  exhibited  statistical  scatter  precluding  a  deterministic  comparison. 
Comparisons  are  made  first  at  the  non-parametric  level  and  then  at  the 
parametric  level  together  with  consequences  in  structural  efficiency  and 
reliability. 

a.       Non-Parametric  Comparison 

Distribution-free  properties  of  the  two  sets  of  data  are 
presented  for  an  objective  comparisons.  However,  because  of  sample  size 
(<25)  is  not  sufficiently  large,  conclusions  should  be  considered  as 
qualitatively  valid. 

The  histograms  of  tensile  strength  for  the  non-preloaded  and 
the  preload  sample  are  presented  in  Figure  13.  Preloading  process  is  seen  to 
have  no  influence  on  the  central  tendency  of  the  strength  data  but  has  a 
preceptable  affect  of  removing  the  lower  weak  tail.  The  reduction  of  the 
lower  tail  is  beneficial  to  both  the  structural  efficiency  and  the  structural 
reliability.  Several  relevant  distribution-free  properties  are  listed  in  Table 
IV. 
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TABLE  IV.  COMPOSITE  STRENGTH  PROPERTIES 

Mean  StdDev        Medium        Mode       Skewness 

No  Preload  46.942  3.07  47.220  47.776        -.27 

Preload  46.418  2.12  46.420  46.424        -.003 

Change  -1%  -31%  -2%  -2%  -100% 

Consistent  with  the  visual  comparison,  column  one  suggests 
no  significant  change  of  the  mean  strength  between  the  two  processes. 
However  column  five  indicates  that  the  nonpreload  strength  skewed  to  the 
weak  tail  whereas  the  skew  is  drastically  reduced  after  preloading.  This  was 
also  consistent  with  visual  identification  of  lower  tail  shift  after  preloading, 
see  Figures  13  and  14. 

b.      Parametric  Comparison 

The  local  load  sharing  model,  [Ref.  3]  provides  the  theoretical 
bases  for  the  Weibull  model  for  the  non-preloaded  samples. 
Characterization  of  the  non-preloaded  strength  data  by  a  two  parameter 
Weibull  model  is  therefore  proper.  Similar  justification  for  the  preload 
sample  is  not  yet  firmly  established;  this  data  can  also  be  fitted  by  a  two 
parameter  Weibull  model  for  expediency  for  comparison,  with  the 
understanding  that  the  procedure  is  ad  hoc.  The  Weibull  parameters 
obtained  through  maximum  likelihood  estimator  for  the  two  sets  of  data  are 
tabulated  in  Table  V. 
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Figure  13.  Histogram  of  Composite  Failure  Load  Results 
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Figure   14.     Weibull  Distribution  Plot  of  Composite  Failure  Load 
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TABLE  V.     COMPOSITE  WEBBULL  PARAMETERS 


shape 

Scale 

Nopreload 

17.48 

48.34 

Preload 

23.45 

47.43 

Change 

+34% 

-2% 

The  Weibull  representations  are  also  consistent  with  the  non- 
parametric  observations,  namely  no  significant  central  tendency  shift  as 
reflected  to  almost  identical  p  values;  but  a  significant  decrease  in  scatter  as 
reflected  in  the  increase  in  a.  The  increase  in  a  affects  structural  efficiency  ( 
at  a  given  reliability  level).  The  magnitude  of  the  potential  improvement  is 
illustrated  in  Table  VI. 

TABLE  VI.    SUMMARY  OF  RELIABILITY  AS  A 
FUNCTION  OF  APPLIED  LOAD  (KGS) 


liability 

Nonpreload 

Preloaded 

Change 

1x10-3 

32.56 

35.33 

8.5% 

1x10-5 

25.02 

29.03 

16% 

1x10-6 

21.93 

26.31 

20% 

For  composite  applications  where  functional  reliability  of 
1x10-3  is  required,  structural  efficiency  goes  up  by  8.5  %.  For  applications 
where  national  security  reliability  is  lxlO-5,  structural  efficiency  goes  up  by 
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16%.  For  composite  applications  where  man-safe  reliability  (  lxlO6  )  is 
required,  structural  reliability  goes  up  by  a  significant  20%. 

Alternately,  an  increase  in  a  can  be  interpreted  as  to  cause  an 
increased  reliability  at  a  given  design  stress  level,  (defined  as  a  fraction  of  the 
mean  strength).  For  example,  The  probability  of  failure  (F)  at  80%  of  the 
mean  decreases  from  1.2%  for  the  nonpreloaded  composite  to  .32%  for  a 
preloaded  composite,  a  delta  of  75% 

Finally,  it  was  pointed  out,  earlier,  that  modeling  the 
preloaded  data  by  the  two  parameter  Weibull  model  was  ad  hoc.  In  fact,  the 
preloaded  data  when  presented  alone  suggest  a  bimodel  trend  (Figure  15). 
Fitting  a  shape  parameter  ,  a  to  the  lower  tail  return  a  value  of  a  =  270  .  This 
large  value  of  a  effectively  provides  for  a  structure  that  would  be  practically 
100  %  reliable  at  any  load  level  below  the  transition  between  the  modes 
(approximately  43  kilograms)  or,  a  100%  increase  in  strength  at  lxlO"6 
reliability.  Thus  potentially,  the  improvement  from  preloading  is  even 
greater.  However,  it  must  be  understood  that  the  number  of  samples  tested 
herein  is  insufficient  to  establish  a  bimodel  system.  This  latter  interpretation 
is  subject  to  confirmation  by  more  extensive  exploration. 
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Figure  15.     BiModel  Trend 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  individual  fiber  testing  it  can  be  concluded  that  Individual 
fiber  Test  system  developed  at  the  Naval  Postgraduate  School,  NPSIT,  is  a 
viable  testing  method  which  produces  accurate  and  reproducible  results. 
The  NPSIT  test,  is  a  one  step  test  system  which  provides  for  the 
determination  of  fiber  diameter  and  ultimate  failure  load.  This  test  method 
minimizes  the  handling  of  the  sample,  thus  providing  a  truer  characterization 
of  the  fiber.  An  additional  advantage  is  that  the  test  was  performed  in  40 
percent  less  time  as  compared  to  using  the  INSTRON  tester  in  conjunction 
with  the  Scanning  Electron  Microscope. 

This  study  clearly  demonstrated  that  preloading  a  composite,  prior  to 
complete  polymerization,  has  a  profound  effect  on  the  lower  tail  of  the 
Weibull  distribution,  hence  the  overall  reliability  of  the  composite.  In 
addition,  the  optimum  preload  level  can  be  independently  mathematically 
estimated,  based  on  fiber  tests  results. 

Further  studies  are  recommended  in  the  areas  of; 

Use  of  the  Micron  Eye  in  determining  fiber  diameter. 

Alternate  composite  gripping  methods. 

The  effects  of  preloading  as  a  function  of  strain  rate,  air 
curing  time  and  different  strains. 

Preloading  composite  strands  at  elevated  temperatures. 

Statistical  modeling  to  account  for  a  multi-model  distribution 
resulting  from  preload. 
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APPENDIX  A.  SINGLE  FIBER  TEST  PROCEDURES 

1.       FIBER  SAMPLE  PREPARATION 

The  individual  fibers  were  mounted  on  cardboard  forms  for  ease  of 
storage  and  handling.  The  samples  were  prepared  in  accordance  with  ASTM 
procedure  D-3379-75  (reapproved  1982)  [Ref  11.]. 

Fabricated  cardboard  mounts.  (3"  X    5/8"). 

Remove  10  inch  strand  of  graphite  form  roll. 

Using  a  magnifier,  carefully  separate  strand 

to  expose  individual  fibers. 

Pick  up  one  end  of  fiber  with  adhesive  tape  and  cut  to 

approximately  twice  the  desired  length. 

Lay  fiber  over  cardboard  mount  and  center  (  Figure  16). 

Affix  the  ends  of  the  fiber  with  tape. 

Place  cardboard  on  storage  rack  and  label. 

Apply  Duco  cement  at  the  ends  of  the  fiber. 


Adhesive 


"0 *F"--' 


3- 


-Q- 


■Gage  Length  =  5  cm — ► 


Figure  16.  Fiber  Mount 
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2.      CALIBRATION  OF  FIBER  DIAMETER  SYSTEM 

Before  the  fiber  diameter  can  be  determined  using  the  diffracted  laser 

beam  associated  with  the  Naval  Postgraduate  School  Integrated  Test,  NPSIT 

a  calibration  curve  must  be  established  . 

The  calibration  curve  used  to  determine  the  fiber  diameter  shows  the 

results  when  comparing  the  fiber  diameter  as  measured  using  the  diffracted 

laser  beam  (NPSIT)  compared  to  the  diameter  recorded  when  using  the 

Scanning  Electron  Microscope  (SEM).  Six  test  specimens  were  mounted  on  a 

15  centimeter  cardboard  mount  and  then    samples  from  each  end  were 

removed  and  carefully  marked.    These  samples  were  then  put  under  the 

Scanning  Electron  Microscope  (SEM)  and  several  diameter  readings  were 

recorded  and  averaged.    The  Scanning  Electron  Microscope  results  were 

used  as  a  standard  for  determining  the  fiber  diameter. 

Figure  17  is  a  graphical  representation  of  the  calibration  result,  with  the 

equation  of  the  straight  line  being. 

Diam  =  .89  *  (NPSIT)  +  .3614  (  1) 

A  linear  curve  fit  proves  most  useful  and  the  results  for  an  exponential 
or  polynomial  curve  fit  are  essentially  the  same. 

It  should  be  noted  that  as  the  fiber  diameter  increases,  the  difference 
between  the  NPSIT  results  and  the  corresponding  SEM  results  also  increases. 
By  reviewing  Figure  6,  this  can  be  predicted  because  as  the  fiber  diameter 
increases  the  spacing  between  the  nodes  (X)  decreases,  and  so  does  the  width 
of  the  associated  dark  areas.  Therefore  it  becomes  much  more  difficult  (and 
less  accurate)  to  find  the  center  of  the  low  intensity  node  using  the 
photocells. 
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Figure  17.  Diameter  Calibration  Curve 


3.       FIBER  DIAMETER  MEASUREMENT  PROCEDURES 

The  fiber  sample  is  first  mounted  in  the  test  by  inserting  between  the 
upper  and  lower  clamps  Then  the  entire  stand  is  adjusted  longitudinally  to 
achieve  the  desire  distance  (L)  from  the  photocells,  (note:  the  spread  of  the 
light  pattern  increases  as  L  increase).  The  photocell  pedestals  are  then 
adjusted  to  an  approximate  position  corresponding  to  the  predicted  distance 
between  the  nodes,  (  ensure  there  is  sufficient  travel  in  the  micrometers  to 
allow  for  left  and  right  movement) .  After  the  pedestals  are  secured  in  place 
measure  the  distance  between  the  photocells  and  record  as  C.  Mount  the 
micrometer  in  place  on  the  photocell  pedestals  and  zero  the  reading.  Now  the 
final  distance  between  the  photocells  can  be  easily  determined  by  adding  C  to 
the  micrometer  reading  Delta  C. 
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The  final  steps  for  determining  the  fiber  diameter  are  performed  with 
the  lab  lights  off,  therefore  complete  familiarization  of  all  test  equipment 
adjustments  is  mandatory.  After  the  laser  is  turned  on  and  the  shutter 
opened,  a  red  light  (beam)  will  appear  (avoid  looking  directly  into  the  laser 
beam).  It  may  be  necessary  to  manually  adjust  the  laser  position  and  the 
transverse  position  of  the  test  stand  to  obtain  the  optimum  diffraction  pattern 
at  the  point  of  interest  on  the  fiber.  Then  manually  adjust  the  height  of  the 
photocells  to  align  the  light  pattern  to  the  center  of  the  photocells.  The 
preliminary  adjustments  are  now  complete. 

To  determine  the  distance  between  the  second  nodes,  each  photocell 
pedestal  is  adjusted  using  the  micrometers.  The  final  adjustments  are  made 
transversely  to  locate  the  position  associated  with  the  highest  resistance 
reading  on  the  multimeter.  Readjust  both  pedestal  separately  several  times  to 
ensure  maximum  accuracy.  After  you  secure  the  laser  power,  record  the 
digital  micrometer  reading  as  Delta  C. 

Calculate  the  fiber  diameter,  a,  from  equations  (2)  and  (3)  of  [Ref  9.]. 

X  =  C  +  DeltaC,  (2) 

a  =  n  (lambda)  L  /  (X/2),  (  3) 

n  =  corresponding  node  (eg  n  =  2), 
L  =  598.6mm.       C  =  115  mm. 
The  accuracy  of  the  fiber  diameter  is  dependent  upon  the  operators 
techniques,  so  five  or  six  readings  should  be  recorded  and  averaged. 
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4.      LOAD  TESTING 

Fiber  failure  load  was  determined  using  either  the  Instron  Universal 
Testing  Instrument  Model  4206  (INSTRON)  or  the  NPSIT  system.  A 
Sensotec  Inc.  50  gram  Load  Cell  Model  31/1435-01  was  used  to  sense  fiber 
load.  The  Machine  Complaince  test,  a  cross  head  speed  of  .1  millimeters  per 
minute  was  best  and  for  the  Fiber  Test  a  speed  of  5  millimeters  per  minute 
(.1  x  gage  length)  was  used. 

INSTRON  testing  is  a  semi-automatic  process  when  used  in  conjunction 
with  the  HP  -  85  computer  and  an  HP  3497A  data  acquisition  system  and  the 
Interactive  Data  Acquisition  software  outlined  in  Appendix  B. 

NPSIT  testing  also  used  the  same  data  acquisition  system  ,  however 
manual  operation  of  the  drive  motor  is  required.  The  Control  Technics 
Corp.  Linear  Actuator  Model  CTC-160  (gear  head  ratio  262:1)  and 
associated  Micrometer  Controller  Model  MMC-10  provide  the  same 
functions  as  the  Instron  cross  head.  The  Speed  Enter  was  set  at  -200  on  the 
thumb  wheel  to  achieve  5  mm/min  and  the  Ratio  Enter  was  set  at  40  for  a 
total  travel  of  45  turns  in  one  minute. 

After  load  cell  calibration  and  system  compliance  tests  are  complete, 
Fiber  testing  was  conducted  as  follows. 

Run  FTST  program. 
Carefully  mount  sample  in  clamps. 
-       Select  DIAM  (NPSIT  only). 

Measure  Fiber  Diameter  (NPSIT  only). 
Select  INPUTS  and  input  initialize  data. 
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Select  ACQUIR,  then  carefully  cut  the  cardboard  mount 
near  the  bottom. 

Remove  slack  in  fiber  (NPSIT  only) 
Hit  program  CONT  and  controller  R/S  simultaneously 
After  fiber  breaks  and  data  acquisition  stops,  stop 
motor  controller. 
Graph  and  store  data. 

With  the  ENSTRON  test  system  fiber  diameter  was  determine  from  the 
NPSIT  first.  Also  this  system  has  automatic  fiber  slack  removal  functions. 


55 


APPENDIX  B.    INTERACTIVE  DATA  ACQUISITION  SOFTWARE 

The  Interactive  Data  Acquisition  Software  (IDAS)  was  designed  for  use 
with  the  HP-85  microcomputer  and  the  HP-3497A  Data  Acquisition  System. 
This  software  consists  of  a  set  of  programs  used  in  determining  the  ultimate 
tensile  strength  and  diameter  of  an  individual  fiber  or  filament.  With  the 
IDAS  it  is  possible  to  create  a  load  cell  calibration  curve  and  to  obtain  output 
in  units  of  grams-force.  A  program  was  also  created  to  determine  the  system 
compliance  in  accordance  with  ASTM  procedures.  Provisions  for  complete 
data  analysis  are  provided  through  data  storage,  retrieval  and  plotting 
routines.  The  software  gathers  data  at  a  rate  of  600  points  per  minute  and 
records  and  plots  Load  vs  Displacement  of  the  fiber  sample. 

The  IDAS  can  be  used  with  the  INSTRON  Model  4206  materials  tester 
or  the  Naval  Postgraduate  School  Integrated  Test  system  (NPSIT).  The 
programs  listed  are  associated  with  the  INSTRON  test  system.  When  using 
the  NPSIT  system,  the  program  steps  that  provide  for  automatic  control  of 
the  INSTRON  (OUTPUT  "704")  must  be  deleted.  Crosshead  control  and 
pretensioning  of  the  fiber  is  accomplished  manually.  A  program  sub  routine 
is  listed  in  the  FTST  module  and  is  used  to  determine  fiber  diameter  using  the 
NPSIT  system.  Functions  keys  are  defined  to  provide  ready  access  to 
commonly  used  sub  routines. 

Mr.  Jim  Nageotte  served  as  the  principle  computer  technician  during  the 
development  of  these  programs. 

Complete  listing  of  the  IDAS  programs  and  a  brief  description  follows. 
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1.  LOAD  CELL  CALIBRATION  MODULE  (LDCALB) 

J^     '       IIM    NAG^TE^nCT    26    36  "LDCALB"   calibrates   the   load   cell   to 

-f,    fPfip  convert  the  voltage  out  into  grams 

4  0    D  I  M    H  <  1  0  0  >  >  E  ■'  1  0  9  •  ..  Z 1  (  1 2  0  >  ■  Z  2 

•v  1  2  Q  "' 
50  DISP  "LORD  '"ELL  CALIBRATION" 
66  CLEAR  ?$'? 

70  DISP  "NUMBER  OF  CALIBRATION 

POINTS    Nl  =  "  Select  tlic  number  of  calibration  points,  (six 

80  INPUT    Ml  is  optimum,  maximum  ten) 

90  FOR    K=l     TO    Kl 

190  DISP     "INPUT    LOAD    LEVEL"  Start  of  loop,  once  through  for  each  data 

110  INPUT    L'K'  point. 

120  W=23    S    E=5000 

1  3 0  R  E  M     :  0 U  f  P  U f    709     ;  "  A R  y  R 1 V T3  Y  D  Input  load  level  applied  to  load  cell. 
5UR0PE2" 

140  ON  TIMER*  1  E  GOTO  240 

150  1=1 

1  6 0  0 IJ  T  P U  T    7  0 ?       "  A  I  0  "  The  3497  A  reads  the  output  voltage  A(I). 
17  0  ENTER    709     =     A  <  I  >  And  excitation  voltage  B(I). 

ISO  WAIT  10 

190  OUTPUT  709   "A! 10" 

200  ENTER  709  ,  B< I  ) 

210  1=1+1 

220  MATT  23 

2 3 0  G 0 TO  16 0 

240  OFF  TIMER*  1 

250  N=I-1 

260  P  1  =  O  i?  P  2  =  0  £  Q  1  =  O  y  Q  2  =  0 

2  70  X1=0  til  X2  =  @  '•!  B0  =  O 
2S0  FOP  1=1  TO  H 

2^0  B0=B6+B ■  I  • 

700  Xl=Xl+fl''  I>'B'  I  ■  .     .  ,.        . 

~  1  0  X  2  =  X 2  ■♦■ '  A  < I  >    B  (  I  )  "'  •'-  2  Fitung  the  data  to  a  linear  function. 

320  P 1 =  P 1 + A  v  I  - 
330  P2=P2+A<  I  V2 

3  4  vj  q  1  =  Q 1  +  B  <  I  > 
350  Q2  =  02  +  B  <  I  >  '-2 
3S0  NEXT  I 

3  70  E?=P1'N 

3S0  D9=P?-N-E9~2 

390  S9=D9'  . 5 

4  00  C  9  =  A  B  S  ••  S  9  -  E  9 1 1 0  @  ) 
4  10  ES=Q1/N 

420  08  =  02  N-E8'x2 

430  88=08^.5 

440  C8= APS  (  <58*'  ES *  1  OO  > 

450  X500=X1'N 

4  60  BO  =  B0''N 

470  U<K>=H'5<K) 

4  80  X6  =  X2-'N 

490  D=X6-J'5<K)A£ 

50O  S 1  (.  K  >  -  D  A  .  5 

510  C3=S1<K>/X5<K) 
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520 
530 
54  0 
550 
560 
5  70 
5  8  0 

5  9  0 

6  0  0 
6  1  0 
620 
6  3  0 
640 
650 
660 

6  7  0 
68@ 
690 

7  0  0 
710 
720 
770 
740 
750 
760 
7  70 
780 

7  9  0 
800 

8  1  0 
820 
8  30 
840 
850 
860 

8  7  0 
880 
890 
900 
910 
920 
930 
940 
950 
960 

9  70 
980 
990 

100 

101 

102 

103 

104 
105 


C3 
C3 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PR 
PF 
PP 
PR 
PP 
PP 
HE 
XI 
X3 
F0 
XI 
Yl 


HE 
XI 


=  1  0  0  *  C  3 
=RBS< C3 


"LORD=" ■ L< 
"MEAN  y=" 

"  C;  f|    ',.'  =  " 

"  r  i.i  i  •  •  i  =  " 


"L  C.i":  UT  PIJT  VOLTRGE 
"  ME  RH  ;..;-"  .;  E9.:  "  'J" 

"  '-,   [j     ',,'  —  "  ;  S'^  ;  M   IJ  II 

"  C  '••'  '  ':>  =  "  ;  C9 
"  EXC  I  TRT I  ON ~ UOLTRGE 
"MERN  l..'="  ;  EQi  "  U" 
"  S  .  D  .  V  =  "  ;  S  8  ;  "  V  " 


"C 
■I  ii 


VI 
Dl 
CI 
C2 
PP 
PR 
PP 
R6 
FO 
G2 
R6 
HE 
PR 
DI 
Dl 
IN 
DI 
IN 
U2 
Ul 
FO 
O  I 
O  I 
O  N 

0  P 
0  P 

0  D 


INT 
INT 
INT- 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
XT  K 

=  0  Of    Y1=0 
= 0  i?  X  2  =  0 
P  K=l  TO  Nl 
=  X  1  +  V  (  K  ) 
=YH-LOO 

=X2+V<K>*2 

XT  K 
=X1/N1 
=  X2^-N1 
=X3/N1 
=  Y1.  HI 
=X2-Xl-2 
=  •'  X  3  -  X 1  *  Y 1  •  -  D  1 
*X1> 


K>;"     ?m 
X5  00  :  " 

S  1  (  K  >  .;  " 

;  C3 


Begins  printing  out  statistical  parameters. 


M 

0 


1 1  ii 


i 


Dl 


Continues  once  for  each  data  point  (load 
level). 


For  least  squares  method,  voltage  =  V(K)  and 
Load  =  L(K) 


INT  "L^Rt'.'  +  B" 
INT  "R=";i:i 
INT  "B=" ;C2 

=  0 

R  K=l  TO  Nl 

=  ^:L(K>-Cl:*y(K>-C2>'v2 

=P6+G2 

XT  K 

INT  "PESIOUAL=" , R6 

SP  "SCREEN  GRAPH  POUT  I  HE" 

SP  "MAX.  LORD=" 

PIJT  Y2 

SP  "MIN.  LORD=" 

PUT  Yl 

=X5< 1 > 

=  X5<  1> 

P  1=1  TO  Nl 

F  X5U>>U2  THEN  U2  =  X5<I> 

F  X5CIXU1  THEN  U1=X5U> 

EXT  I 

PINT  " X  MPX  •;  1  000V/ V )  =  "  ;  U2 

1000 

PINT  " X  M I N  < 1O0O  V / V >  =  " i U 1 

1  O  O  0 

ISP  "INPUT  MAX  FOR  X=" 


After  all  load  levels  are  recorded,  the  program 
calculates  and  returns  the  best  fit 
coefficients,  A  and  B. 


Determines  the  goodness  of  fit  which  is 
expressed  by  the  "residual". 


The  graph  routine  plots  the  load  in  grams 
versus  voltage  out  of  the  load  cell.  This 
graph  appears  on  the  computer  screen  and  is 
printed  by  the  thermal  printer. 
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066  INPUT  U2 

070  DISP  "INPUT  MIH  FOR  X=" 

088  INPUT  U 1 

090  FOR  1=1  tq  Nl 

1  00  B  1  •  I  >  =X5  <  I  '>  1  1  000 

110  B2< I >  =S1  '-  I  .  *  1  000 

120  NEXT  I 

1  30  91  =  160  '<:U2-U1  ) 

140    S2=140/<;Y2-Y1> 

159  FOP  1=1  TQ  5 

160  GRAPH  <k    L.CLEAP 
170  LDIP  0 

1  SO  SCALE  -48  . 208 . -36 ,156 
190  X fl  X  I S  0 , 3  2  0 ,  1 6  0 
200  VOX  IS  0.280.1 40 
210  PEN  1 
220  PENUP 

230  Wl=160  'OJ2-U1 > 
240  W2=146^<Y2-Y1  '■' 
250  FOP  1=0  TO  5 
260  U  =  U  1  +  (  U2-U  1  f  "'5*  I 
270  P3=- U-Ul )*H1 
280  MOVE  P3, -  1  @ 
290  LABEL  VflL$<U> 
300  NEXT  I 
310  FOR  1=0  TO  5 
320  Y=Y1+<Y2-Y1  >'5*I 
330  Q3=<Y-Y1 >*te2 
340  MOVE  -23-Q3 
350  LABEL  VflL*<Y> 
360  NEXT  I 
370  MOVE  40,  -2-: 
380  LABEL  "  1  OOOtV.'ou  t  /Vexc  " 
390  DEC 

400  MOVE  -29, 25 
410  LDIP  90 
420  LABEL  "F0RCE<<3m>  " 
4  30  LDIP  0 
440  POP  1=1  TO  Nl 
450  Zl< I >  =  < B 1 < I >-Ul > *W1 
460  Z2<I>=<L<I>-Yl)*W2-2 
4  70  MOVE  Zl(I)iZ2(I) 
480  LABEL  "*" 
490  NEXT  I 
5O0  G1=C1  1OO0 
510  FOR  I=0  TO  lOl 
520  U  =  U  1  +  <  U2-U  1  >  ■'  1  @0* I 
530  Y=G1*U+C2 

540  IF  Y<Y1  OR  Y>Y2  THEN  1580 
550  Zl  <I>  =  <U-UJ  )U\\ 
560  Z2<I>=<Y-Yi>*M2 
570  PLOT  ZKI),  Z2-:i) 

580    NEXT     I  This  routine  plots  the  graph  on  the  HP 

590    COPY  Plotter. 

6O0    DISP     "DO    YOU    WANT    A    HARD    CO 
PY    ON    PLOTTERS <Y/N>" 
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1610 
1620 
1630 
1640 
1650 

1660 

1670 


168  0 

IF.  98 

17O0 

1710 

1720 

1730 

1  740 

1750 

1760 

1778 

1780 

1790 

1800 

1810 

1820 

1S30 

1840 

1850 

I860 

1  870 

1880 

1  850 


1  9  O  0 

1910 
1920 
1930 
1940 
1950 
I960 


19  70 

1930 
1990 

2000 
2010 

2820 
2830 
2848 
2850 


INPUT  K8$ 

IF  K8t="H"  THEN  2176 

PRINTER  IS  10 

CONTROL  10=5  i     48 

OUTPUT  10   "I N  > SP 1 ,  I P24O0 ,  ! 

608..  8808  ■  6900  i  " 

OUTPUT  10   "  SCO  ,  1  00O  •  0  .-  1  000 

i. 

OUTPUT  I  0   " PUO • OPDO •  1 800  >   1 

O0O  •  1  800  ,   1  COO  •  O  •  0  .•  OPU  " 

W  3  =  1  0  0  0 ■-'  <  U  2  -  U  1  > 

W4=lPifu3' ■<  Y2-Y1  > 

OUTPUT  10   "SIO  2, O  3 i TL1 .5 

.  0" 

FOR  1=8  TO  5 

V  =  Y1  +  '.:Y2-Y1  >/5*I 

Y4=0'-Y1 >*W4 

Y'=INT<  Y> 

Y4=INT<Y4> 

0 U T P U T  1 O  ;   " PR  8 ,  " , V 4 ,  " Y T ,  " 

OUTPUT  18  i"  CP-5  ,  -O  .  07  .;  LB  "  , 

Y;CHR*<3> 

ne::t  I 

FOR  I=0  TO  =■ 
U  =  U  1  +  <  U  2  -  U  1  >  *  I  y  5 
U4='U-U1 >*W3 
U4=INTMI4> 
U=INT<U) 

OUTPUT  10  i  "  PR  "  ,  U4  ,  "  ,  O .;  XT;  " 
OUTPUT  10  i "CP-1 .3,-1 iLB" jU 
, CHR*<3> 
NEXT  I 

OUTPUT  10  .;  "SI.  30,.  42" 
OUTPUT  10  .;  "FR400,0.;CF-2,-2 
.  3;LB1800*Vout.'V<?xc"  ;CHR*<3 
> 

OUTPUT  10  .;  "PR0,460;D10,  1;C 

P-2.6,2.6;LB  FORCE  ( <?rn  >     ";C 

HR$<3) 

OUTPUT  10  ;"OI,PU" 

FOP  1=1  TO  HI 

Z1<I)=<B1<I)-U1)*M3 

Z2< I>=<L<I>-Y1 >*W4 

Zl<  P  =  IHT(Z1(I.)) 

Z2<  I>  =  INT<Z2< I >> 

OUTPUT  1  O  :.  "  PU  ,  "  ;  "  PR  "  ,  Z  1  <  I  > 

, Z2< I> ; "CP-0. 1 , -8 .2; LB*" ; CH 

R*<3> 

NEXT  I 

OUTPUT  10  ; "PU" 

FOR  1=0  TO  101 

U  =  U 1  +  < U2-U 1 > / 1 00* I 

Y=G1*U+C2 

IF  Y^Yl  OR  Y>Y2  THEN  2060 

ZKI)  =  ai-Ul  >*W3 

Z2CI >=<Y-Y1>*W4 

OUTPUT  10  ; "PR";Z1<I >,Z2< I) 

;"PD" 


60 


2060  NEXT  I 

2070  OUTPUT  10  .;  "PUO,  900..  100  •  900 

if 

2 680    OUTPUT     10     .;  "PU"  The  user  can  enter  a  label  or  legend   on  the 

2O90    OUTPUT     10     ;  "SI.  22,     38"  finishedplot 

21O0    DISP    "ENTER    THE    LEGEND.     ENT 

ER  'O'  TO  EXIT" 
2110  INPUT  P7* 
2120  IF  P7*="U"  THEN  2150 
2130  OUTPUT  10  .;  "CP;LB"  ;P7$  =  CHP* 

2140     GOTO     2100  ■•inrAinVnHo 

2150    PRINTER     IS    2  LDCALB  ends. 

2160  0I3P  "  END  LDCHLB" 

2170  END 
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2.  MACHINE  SYSTEM  COMPLIANCE  MODULE  (CTST) 

18  !  CTST37 

20  '        JIM  NAGEOTTE          JRN     13    87      The  system  Compliance  module  or  "CTST" 

■-■y  •_  measures  the  compliance  or  displacemcni  of 

4 0  L L E  AR  the  system  without  a  fiber  sample. 

59  disp  @  disp  a  cusp  @  disp  e 

DISP  "  CT8T87" 

60  DISP  @  DISP  "   USE  THE  FUNCT 

ION  KEYS . " 

70  PRINTER  IS  2  tu-  j  ■       r 

L:     r-t.Vf^^cr—     ,.,_..     ._..,__.     _.   ,   This  program  is  menu  driven  using  function 

80    DIN    P£45o>^45y,.,CK450;...  Zl<   kcysVacccss  the  sub  routinc£  should 

90    PRINT^IF^FROGRAM    HALTS,     TYP  f*c.culi°n  ^„  relurn  t0  the  mcnu  by 
E.'CGNT     108'     TO    RETURN    TO    ME    ypmg  " 

NU    " 

180  PRINT 

1 10  GOTO  130 

128  CLEAR 

130  ON  KEY#  •*, "END"  GOTO  220      n  _  JU 

140    ON    KEY*     1,  "INPUTS"    G0SU6    2 4 0  The keys are dcfined nere- 

150    0 N    K E  Y #    2  ,  "  fl C Q U  IP."    G 0 S U B    53 O 

160  ON  KEY#  3," DISP"  GOSUB  1118 

170  ON  KEY*  6,"PD_DAT"  GOSUB  156 

0 

ISO    ON    KEY*    7     "FLOTTP"     GOSUB     136 

0 

190  ON  KEY#  5.  "SA'..'_DAT"  GOSUB  17 

10 

200  KEY  LABEL 

210  GOTO  130 

220    CLEAR    i?    DISP    "END    CTST87"    5     The  program  will  continue  to  loop  on  the 

END  menu  until  a  function  key  is  depressed. 

2  30     ! 

240     !        **     INPUTS     **■ 
250    CLEAR 
260     !     D I  SP  "  ENTER    DATA    RRTECpts'm  The  "INPUTS"  sub  routine  is  required  first. 

i  n  )  "  Input  the  variables  and  information  needed 

£78     !  I NPUT    P  along  with  the  values  of  A  and  B  obtained 

280    R  =  6O0     !     PRESET    DATA    RATE  from  "L  DC  ALB". 

290    DISP    "ENTER    ELAPSED    TIME    (3D 

OP  45:'" 
3Pih  INPUT  EO 
310  CLEAR  @  DISF  "ENTER  SPEED  IN 

NM'MIN" 
320  INPUT  S 

325  PRINT  "SPEED  =  "  .;  S  i  "MM.-M  I  N" 
33«  E1=E0*1O00 

340  DISP  "ENTER  DATA  FILE  NAME" 
350  INPUT  N* 
360  CLEAR 

376  DISP  "P  =  fi*:V  +  B"  •?  DISP 
388  DISP  "ENTER  'A'  FROM  LDCflLB 

LISTING" 
382  INPUT  C3 
480  DISP  @  DISP  "ENTER  'B'  FROM 

LDCflLB" 
410  INPUT  B 
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420OLERP709 

4  7Q  Hi  UP!  IT    7flQ       "Hi  10"  Tnc  IIP  3497  A  records  the     excitation 

446  ENTER    709.;     W  voltage  applied  to  the  load  cell. 

450  DISP    "EXCITATION    VOLTAGE    I  S "  The  equation  that  translates  the  load  cell 
;  im,  voltage  into  physical   units  of  grams   is 

4  g  g  h  =  R  G  S '  f  7    U  ')  adjusted  for  the  present  excitation  voltage. 

4  70  R  =  I  N  T  < H * 1  0  0  0  >  -  1  0  0  0 

480  DISP     "  fi=  "     H  "INPUTS"  is  terminated  by  returning  to  the 

4?0  DISP     "B="B  menu. 

5  0  0  C  0  P  V 
505  CLERP 
510  RETURN 

520  !  The  data  acquisition  sub  routine  is  title 

5  3  0  '        *  *    H C  0  U  I  F        t  i  " ACQUIR".  Acquir  checks  die  time  duration 

540  IF    E0>8    THEN    598  selected  and  if  it  finds  it  to  be  not  greater 

550  BEEP    £    BEEP    @    BEEP    (?    CLERR        than  zero,  prompts  the  operator  to  run 

560  DISP    "'INPUTS'     MUST    BE    PUN    F  "INPUTS"  and  returns  to  the  menu. 

IRSTi  " 

570  RETURN 

598  CLERP 

600  BEEP 

610  FOR    1  =  1     TO    ,     i?    DISP    i*    NEXT    I  The  program  pauses  and  waits  for  the 

626  D  I  SP    "  CHECK     IEEE    BUTTON  •     MQU  operator  to  hit  "cont".  At  that  time  a  timer 

NT    SAMPLE    AND    BURN    "    @    DISP     is  started  and  data  is  acquired  at  a  rate  of  600 

630  DISP    "HIT    CONT    WHEN    READY"        data  points  per  minute. 

640  PAUSE 

650  CLERR    &    DISF 

653  DISP    "FINDING    ZERO    LORD" 

655  OUTPUT    704    ;HK13,".:S 

658  OUTPUT    709     :  "AP.VR1  VT3VD5VA8A 

666  OUTPUT    789     ,-flIB"  Notc:      OUTPUT    "704"    provides    for 

---  --'-'--.  -  automatic    control     functions    to    the 

|~*  p"a+L  +  R  '  IN STRON  tester.    These  functions  provide 

|g3  |  F    p  < 0    T HEN    GOTO    678  for  a  slight  tcnsion  .in  the  fiber  pr'°r  l°  lh° 

664  OUTPUT    7m4     ,  "P"  start  of  data  acquisition. 

665  WRIT    588  .    ,  r- 

666  OUTPUT    784     ."Km"  Thc  starl  of  lhe  timer  ,s  dcfincd  at  mc  zcr0 

667  GOTO  668  load. 
678  CLERR  @  BEEF  @  DISP  G  DISP  " 

MOVING  UF" 

675  0UTPUT  704   "K3" 

680  OUTPUT  709  ; "AI0" 

685  ENTER  70S"  .-     L 

687  P=A*L+B 

688  IF  P<.06  THEN  GOTO  680 
691  OUTPUT  704  i  "K21" 

694  BEEP  S  DISP  "TIMER  ON" 

695  ON  TIMER*  1,E1  GOTO  768 
788  1=1 

718  OUTPUT  789  .:  "A  19" 

728  ENTER  709  ;    L(I) 

736  1=1+1 

740  WRIT  23 

750  GOTO  710 

760  OFF  TIMER*  1 

765  OUTPUT  704  ; "KO" 

770  1=1-1 
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789  WRIT  50 

790  CLERR 

899  FOR  T=l  TO  7  @  DISP  @  NEXT  T 

810  BEEP 

8  30  DISP 

840  CLEAR  @  DISF  "  WRIT" 

850  W=I 

860  PRINT  \k    PRINT  N;"DRTR  POINTS 

M 

865  PRINT 

8 6  6  0  U T  P U T  70  4   " K 1 " 

870  FOR  1=1  TO  N 

a 89  P'.  I  >  =  A$L<.  I  >+B  The  voltage  is  converted  to  grams  force 

8?o  next  i  m)l 

910  P  2  =  P  < 1 > 

320  P1=P<  1  '> 

930  FOR  1=2  TO  h 

940  IF  PCI)  >  P2  T  HEN  P2  =  F  (.  I  > 

950  IF  P  ' I  >  <P  1  T  HEN  PI  =P  <  I  > 

970  NEXT  I 

980  PRINT  «'p' MR:0  =  "  .;  P2 

990  PRINT  "P<MIN>=" :P1 

1000  PRINT 

10S0  CLERR 

1090    RETURN  "ACQUIR"  terminates  by  returning  to  the 

1  IQm  !  menu. 

1110  !    **  DISP  GRRPH  ON  CRT  ** 

1120  PRINT  "  "  ,Nt  5 

PRINT 

1130  CLERR  l?  DISP  "INPUT  P<MflX> 

FOR    Y    R X  I  S  '  The  function  key  denned  as  "DISP"  calls  the 

1  140  INPUT    P2  sub  routine  which  plots  the  data  on  the 

1150  P  1  =  0  screen  and  on  the  thermal  printer. 

1160  GCLERR 

1170  GRRPH 

1 180  LDIP    0 

1199  SCALE  -48,  298  ■  -36..  156 

1 200  XRX  IS  0 , 32 .  0 •  160 
1210  YRXIS  9,28, 9- 140 
1220  PEN  1 

1230  PENUP 

1  240  S  1  =  160-'  <.  EO .'■  60  >  /S 

1250  S2=140'P2 

1260  FOR  T=l  TO  5 

1270  X=S*T/10 

1 280  X= I  NT  < X* 1 900  > / 1 OOO 

1290  X2=X*S1 

1300  MOVE  X2..-10 

1310  LABEL  VAL*<X> 

1320  NEXT  T 

1330  FOR  T=0  TO  5 

1340  P  =  P2--5*T 

1350  P3=P*S2 

1360  MOVE  -23. P3 

1370  LRBEL  VflL*<P> 

1380  NEXT  T 

1390  MOVE  49,-23 

140O  LRBEL  "DISP  (mui)" 
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14  10 

DEC 

1420 

MOVE  -29,25 

1430 

LDIP  90 

1  4  40 

LRBEL  "FORCE  •'?(»>" 

1450 

LDIP  0 

14  60 

FOR  1=1  TO  N 

147  0 

o<  i>=s*<:i/p> 

14  80 

Zl< I >=D< I >*S1 

1  4  3  0 

Z2U  >=p<:  I  >  tS2 

1  5  0  0 

PLOT  Zl<  I  >  ■ Z2<  I  > 

1510 

NEXT  I 

1520 

COPY 

1  5  3  0 

CLEAR 

1540 

RETURN 

1550 

1 

1  5  6  0 

!     ft   PEhD  DATA   ** 

1570 

CLEAR 

1580 

DISP  "IH8EFT  DATA  TAPE" 

1590 

DISP  "HIT  CUNT" 

1  6  0  0 

PAUSE 

1  6  i  0 

ASSIGN*  1  TO  Nf 

1620 

READ*  1  !  N,E0,S 

163  0 

FOR  1=1  TO  N 

1640 

READ*  1  i  FU' 

1650 

NEXT  I 

1  660 

ASSIGN*  1  TO  * 

1670 

DISP  "  DATA  READ" 

1630 

R  =  6  O  0 

1690 

BEEP 

1700 

RETURN 

1710 

'    **   SA"_OAT   ** 

1720 

CLEAR 

1  730 

IF  F "  <  3  ">  <  >  O  THEN  GOTO  179 

1  740 

BEEP  @  BEEP  @  CLEAR 

1750 

DISP  "NO  DATA  IN  MEMORY" 

1760 

DISP  "HIT  CON! " 

1  7  70 

PAUSE 

1780 

RETURN 

1  790 

DISP  "INSERT  DATA  TAPE" 

1  8  0  0 

DISP 

1810 

DISP  "HIT  CONT" 

1820 

PAUSE 

1630 

CLEAR  @  DISP  "WAIT" 

1840 

M=<N+1 >t8+50 

1350 

H=INT<M/256>+l 

I860 

CREATE  N*,H,256 

1870 

ASSIGN*  1  TO  Nf 

1830 

PRINT*  1  ;  N,E0,S 

1890 

FOR  1=1  TO  N 

1900 

PRINT*  1  ;  PCI) 

1910 

NEXT  I 

1926 

ASSIGN*  1  TO  * 

1930 

DISP  "DATA  SAVED" 

1940 

RETURN 

1950 

i 

D(I)  is  the  displacement 


"DISP"  terminates     by  returning  to  the 
menu. 


The  load  data  is  read  in. 

The  number  of  data  points  saved,  the  time 
the  test  ran  and  the  speed  are  inputted.. 


A  prompt  and  beep  indicate  the  data  has  been 
read  in. 

The  sub  routine  is  terminated. 

The  "SAV-DAT"  routine  is  just  the  counter 
part  of  the  above  routine.  It  first  checks  to 
see  that  there  is  data  in  the  memory  to  be 
saved. 


The  memory  space  is  computed  and  the  file 
opened.  The  data  is  then  saved.  It  is  in  the 
same  form  as  that  read  in  by  "RD-DAT". 


The  sub  routine  is  terminated 
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I960     '  **    PLOTTER    **  _     „_.  r,^rnD„  r      ,.      .    .    . 

1Q70    i'LEhP  PLOTTER    function  is  design  to  create 

1930    D  ISP    "ENTER    Y    MRX    FOR    Y    fl  x  j  *  P^t  on  the  HP  Plotter. 

1999  INPUT  Y2 

2000  DISP  "ENTER  Y  MIN  FOR  V  FIX  I 

c  it  Select  the  maximum  and  minimum  values 

2  010     I  N  P  U  T    Y  1  f°r  ^c  X  and  Y  axis  on  the  plot. 

2020    DISP    "ENTER    X    MAX    FOR    X    HXI 

i~-  it 

2030    INPUT    X2 

2040    DISP     "ENTER    X    MIN    FOR    X    RXI 

5  "  This  plot  will  be  on  the  HP  plotter  similar 
2  0  50     INPUT    XI  to  the  plot  on  the  CRT  screen. 

2  0  6  0    S  3  =  1  0  O  0  /  (.  X  2  -  X  1  > 

2  0  70    S 4  = 1 0  0  0 / <  Y £ - Y 1 > 

2O80  PRINTER  IS  10 

2090  CONTROL  19,5  :  43 

2 1  00  OUTPUT  lO  .;  "  IN .;  SP 1  ;  I P2400 ,  1 

6  O  O • 8 3  O  O ■ 6900 :  " 

2110    OUTPUT     10     ; "SCO, 1000,0, 1009 

n 

2  1  20    OUTPUT     lO     •  " PUO , 0PDO ■  1 909 ,  1 

000 , 1 000 i 1 00O      0,0, OPU " 

2130  OUTPUT  lO  : "SI9  2,0. 3; TL1  5 

,  0  " 
2140  FOR  1=0  TO  5 
2150  Y=INT^  Y1  +  02-Y1  >.'5*I> 

2 1 69  Y4= I  NT  C  <  Y-V 1 > *S4  > 

2170  OUTPUT  10  .;  "  PR  0  ,  "  ,  Y4  ,  "  YT  ;  " 
2130  OUTPUT  10   "CP-5, -O . 07; LB" ; 

Y;CHR$<3> 
2190  NEXT  I 
220O  FOR  1=1  TO  5 
2210  X  =  X1-KX2-X1  >/5*I 
2220  X4=INTOX-Xl >*S3> 
2230  OUTPUT  1  O  }  "  PR  "  ;  X4  ,  "  ,  0  .;  XT  ;  " 
2240  OUTPUT  10  . "CP- 1 . 3 , - 1 ; LB" .: X 

2250    NEXT    I 

2  2  6  O    0  U  T  P  U  T     10    .;  "  S  I  .  3  O ,  .  4  2  "  The  lables  for  the  axis  are  put  in  place. 

2279  OUTPUT  10  •■  "PR220,O  .CP4, -2  . 

3 j LBDISPLRCEMENT  Cmm> " ; CHR* 

<  3  • 
2230  OUTPUT  19  ;  "  PRO  ■  469  .;  D  19,1  .;  C 

P-2  6,  2.  69;  LB  FORCE  <gm>"  .;  C 

HR*' 3> 

^t«0  p,^p     t  =  i     TO    N  Here  displacement  is  represented  by  D(I) 

2310  D<n=StI    P  while P(I)  is  load. 

2320  Zl<  I  >  =  INK  <  D<  I  )  -X  1  >  *S3  > 

2339  Z2< I>  =  INT<  • .  P< I >-Yl )*S4> 

2340  OUTPUT  19  . "PR"  ,21<I >,Z2<I ) 
.  „pD„ 

2359  NEXT  I 

2369  OUTPUT  19  , "PU9 , 999, 1 90 , 900 
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2379  OUTPUT  19  : "PU" 

2388  OUTPUT  19  ;"?.j     2"-,  70- 

2390  DISP  "LOCATE" PEN^ FOR  LEGEND 

"  After  the  plot  is  completed,  the  pen  moves 

2400    DISP    1?    D  I  S  F     "ENTER    THE    L  E  G  E t0  l^c  uPPcr  'e^  corner  where  a  legend  may 

NO    YOU    WflNl     "  be  created. 

2410    DISP     "ENTER    A     '0'     WHEN    DONE 

H 

2420  INPUT    P* 

2430  IP  P$="0"  THEN  24*fl 

2440  OUTPUT  10  , "CP; LB" :P$; CHR$( 

3  > 

d 4 5 0  ij 0 TO    2  4  0 0  Thc  program  is  completed  and  returns  to  the 

2460  PRINTER     IS    2  menu. 

24  70  RETURN 

2535  CLEAR    709 
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3.  QUADRATIC  FIT  MODULE  (QFIT) 

-  QFIT    fits  a  quadratic  equation  to  the 

--.  fl  j  |  j  ^    ^  h  r  f  n  j  T  t  compliance  data  obtained  from  "MCTST". 

30  ! 

46  !   FIT  the  MCTST  DATA  BY  fl  Q 

UADPATIC  equation 

50  ! 

60  CLEAR 

70  PRINTER  IS  2 

80  SHORT  P< 459 >  ..  Z 1  < 450 >  ,  Z2<  459  ) 

90  DIM  E < 2 , 2 > , F <2 ) , G< 2 > , D (456 > 

100  BEEP  <k    DISP  "INSERT  DATA  TAP 
E  " 

110  DISP  "INPUT  FILE  NAME" 

120  INPUT    \\t  Input  file  name.  Use  the  file  stored  by 

130  ASSIGN*  1  TO  N*  "MCTST". 

140  READ*  1  ,  N-E0.S 

150  FOR  1=1  TO  N 

160  READ*  1  ,  P<I) 

170  NEXT  I 

ISO  ASSIGN*  1  TO  * 

190  P=60O  !  DATA  SAMPLING  RATE<p 

n  t  s  ■■'  m  i  n  > 

200  FOP  1=1  TO  N 

210  0'  I  >  =  ItS-'P 

220  NEXT  I 

230  DISP  "DO  YOU  WANT  TO  DISP  TH 

E    GRAPH    ON    SCREEN    (Y/N)?"  A  graph  or  plot  at  this  point  would  give 

24  0  BEEP  results  identical  to  the  output  of  "MCTST". 

250  INPUT    fl$ 

260  IF  fl$="Y"  THEN  GOSUB  1020 

270  ' 

280  DISP  "DO  YOU  WANT  A  HARD  COP 

Y  ON  PLGTTEP?<Y''N>" 

290  INPUT  A* 

390  IF  A*="Y"  THEN  GOSUB  1600 

310  ! 

320  DISP  "0  FIT  SUBROUTINE" 

330  DISP  "COMPUTING   ABOUT  2  MIH 

i  I  j  E  S    " 

7  4rt  E'Q.0>=N  The  "QFIT"  sub  routine  requires  about  two 

7c1.  c  .'  n    1  \-r*  minutes.  The  exact  time  depends  on  the 

360  E  <  9    2  >  =9  number  of  data  points  to  be  iittcd. 

3  70  EU.1)=0 

380  E<1,2>=0 

390  E-'2,2)=0 

400  FOR    1=21     TO    N  Initializing  the  matrix 

4 1 O  E  < 0 , 2  > =E  <  O , 2 >  +P < I > -2 

420  NEXT  I 

430  FOR  1=1  TO  H 

440  E  ■-.  O  ,  1  >  =E  <  0  ,  1  >  +P  <  I  ) 

441  E  <  1  *  2>=E  (.  1  /  2  )  +P  <  I  >  A2*  <  P  <  I  >  *N  Least  square  method 
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520  E'  2 ,  2  >=E<2,  2>+P<  I  )*2*<P<  I 
530  NEXT  I 


535    Ed  •  ! 
54  0    E  ■:  1  ■  2 


=  E  • '  0  •  2  > 
=  E<1 .2> - N 
550    E<-:2,2''=E<2,2;'    N 

5  7  0    G  <  0  >  =0    i?    G  ••  1  >  =  A    i*    G  '•  2  ■■  =  0 
b90    FOR     1  =  1     TO    H 

b  0  0    G  (  9  >  =  ij  (  0  >  +  0  ■  I  ) 

630    G  •:  1  '  =G  (  1  >  +Dd  )  * P  •  I  > 

6  6  O    ij  <  2  >  =  G  <  2  >  +  D '  I  ■•  *  P  ••  I  )  's  2 
6  70    NEXT    I 

680  DISP  "MRTPi:  OPERATION" 
6  9  0  E  <  1  ■  0  >  =  E  •  0  j  1  > 

7 00  e  <:  2 ,  o  >  =e •: 0 . 2  ■  By symmetry 

710  E<2,  1>=E<L2> 

720  MRT    DISP    E 

730  MAT    F  =  SYS<E,G;> 

740  PRINTER  IS  2 

750  PRINT  "        "  ;N|-  e  PRINT 

760  PRINT     "P  =  E  +  F*0ELTR  +  G*DELTR^2Tnc  bcsl  fit  coefficients  E,  F,  and  G  are 

"  returned  and  printed  out. 

770  PRINT     "E="  .:F>:0) 

780  PRINT     "F=" ;F<1> 

799  PRINT    "G=" : F<2) 

800  DISP  "P=E+F*DELTA+G*OELTA^2" 
810  DISP  "E  =  "  =  F' 0) 

820  DISP  "F  =  "  :F-  1> 

8  30  DISP  "G=" :F<2> 
S40  MRT  PRINT  F 
850  E=F<0> 

860  F  =  F-:l) 

870  G  =  F<:2) 

880  CLEAR  @  DISP  "COMPUTING" 

890  FOP  1=1  TO  M 

9  0 O  D  v  I  >  -E +F  *  P  (  I  >  +G*P  <  I  >  " 2  The  coefficients  are  then  used  to  calculate  the 
9  1  y  N  E . ,  T  estimated  curve  of  displacement  as  a  function 
?20  BEEP  of  load  PG). 

930  DISP  "DO  YOU  WANT  R  COPY  OF 
ESTIMATED  CURVE  ON  THE  SCREE 
N  ?<Y^N>" 

lit    TFFfiI-»v«     TMPW    m-iiP    m->n  A  plot  of  this  curve  ,  shows  a  smooth 

960  DISP%0  YOlfilRNT^H^^ 

Y  OF  ESTIMATED  CURVE  ON  PLOT 

TER-'-Y-N  V 
970  INPUT  A* 

980  IF  R*="Y"  THEN  GOSUB  1600 
990  BEEP 

lOOG    DISP    "END    OF    QFIT"  n 

1010    END  The  main  program  ends. 

1020  !   **  SCREFN  PLOT  ** 

1030  DISP  "DISPLAY  THE  RESULT  ON 

THE    ^ PREEN"  ^nc  screen  plot  routine  is  the  same  as  that 

1  0 4 0    P P  I  N  T  "  "  "  i  N$  uscd  m  "MCTST"  and  plots  the  graph  on  the 

1  0  5  0    P  R  I  N  T  CRT  screen  and  thermal  printer. 

1060    Y2=P(1> 
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0  70    V 1=P< 1  ':« 

080    FOR     1=2     TO    N 

030     IF    P<I>>Y2    THEN    Y2  =  P<n 

100     IF    PdXYl     THEN    Yl=P<n 

110  NEXT  I 

120  PRINT  "V  MPX=" ;Y2 

130  DISP  "INPUT  Y  MRX  FOR  Y  RXI 

140  INPUT  Y2 

150  PRINT  "Y  MIN=";Y1 

160  DISP  " INFUT  V  MIN  FOR  Y  RX I 

170  INPUT  VI 

130  DISP  "ENTER  X  MRX  FOR  X  RXI 

C  " 

130  INPUT  X2 

200  DISP  "ENTER  X  MIN  FOR  X  RXI 

2io  Input  xi 

220  GRRFH 
GCLERR 
LDIR  0 

SCALE  -48..  203  ,  -36,  156 
XRXIS  9,32,0, 166 
YRXIS  O, 28, O ■ 140 
PEN  1 
FENUP 
X1=0 

31  =  160  ■•''  (  X  2  -  X  1  • 
S2=140/<Y2-Y1> 
FOR  1=1  TO  5 
X  =  X2-'5*I 
X3=rX-Xl  >*S1 
MOVE  X3,-10 
LRBEL  VflL*<X) 
NEXT  I 

FOR  I=0  TO  5 
Y=Y1+<Y2-Y1 >/5*I 
Y3=<Y-Y1>*S2 
MOVE  -23, Y 3 
LRBEL  VfiL*<Y> 
NEXT  I 
MOVE  40 ■ -23 

LRBEL  "DISPLACEMENT  (iBm)" 
DEC 

MOVE  -29,25 
LDIP  90 

LRBEL  "FORCE < ?m>" 
LDIR  O 

FOR  1=1  TO  N 
ZHI)  =  (D(I  )-Xl  >:*S1 
Z2CI )=(P(I >-Yl >*S2 
PLOT  ZKI><Z2<I) 


2  30 

240 

250 

260 

270 

2  8  O 

230 

3  0  0 

310 

320 

330 

34  0 

3  5  O 

360 

370 

380 

33  O 

400 

410 

420 

4  30 

440 

450 

460 

470 

4S0 

430 

500 

510 

520 

5  3  O 

540 

550 

560 

5  7  0 

580 

530 

60O 

The  X  axis  or  displacement  is  represented  by 
D(I)  while  the  Y  axis  or  load  is  represented 


NEXT     I  bypW- 

COPY 

RETURN  The  plotter  routine  plots  the  graph  on  the 

i  HP  Plotter. 

!  **  FLUTTER       ** 
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610  UISP  "ENTER  Y  MAX"  <*     INPUT 

Y2 
620  DISP  "ENTER  Y  M I N "  S  INPUT 

VI 
636  DISP  "ENTER  X  MRX"  @  INPUT 

640  DISP  "ENTER  X  MIN"  5  INPUT 

XI 
6  5  0  S  3  =  1  9  0  y  '■'  <  X  2  -  X  1  > 
•5  6  0  S  4  =  1  0  0  0  /  (  Y  2  -  Y  1  '•' 
670  PRINTER  IS  10 
680  CONTROL  10,5  .:  48 
690  OUTPUT  10   "  IN; SP1 ;  IE2400,  1 

6  O  O  ,  8  8  O  O  •  6  9  O  0  .;  " 
70O  OUTPUT  10  , "SCO, 100O • 0, 1000 

710  OUTPUT  10   "PU0,0RD0, 1000, 1 

0  0  O  t    1  0  0  O  •  1  O  0  O  •  O  •  0  ,  O  P  IJ  " 

720  OUTPUT  10  •  "SI0.2,0.3;TL1 .5 

>   O" 

770  POP  I=0  TO  5 

7ay  Y=INT-  Yl  +  (  Y2-Y1  ':  -5*1  > 

750  Y4=INT< v Y-Yl >*S4> 

760  OUTPUT  10  :"PR  0, " , Y4 . "YT; " 

770  OUTPUT  10   "CP-5,-0  07:.I_B"; 

Y  ;  CHRf  <!3> 
7SO  NEXT  I 
790  POP  1=1  TO  5 
8O0  X  =  X1  *■<  X2-X1  ,'  ''5*1 
810  X4=<X-X1>*S3 
820  X4=INT«'X4> 

830  OUTPUT  lO  ."PR  "  ,  X4  ,  "  ,  O  .;  XT  ;  " 
840  OUTPUT  lO  , "CP-1  3,-l;LB";X 

;CHR$<3> 
850  NEXT  I 

860  OUTPUT  10   "SI. SO,. 42" 
870  OUTPUT  lO  ;  "PR370,0.;CP-2, -2 

.  3, LBDISPLRCEMENT  <mm>";CHP 

*  (.  3  > 
880  OUTPUT  10  ■"  PRO  ,  460  ;  D  I  0  ,  1  :  C 

P-2.6,2.68;LB  PORCE  (qni)".;C 

hp$-:3) 
890  output  10  = "di ;pu" 

900  FOP  1=1  TO  N 

910  Zl<  I  >  =  I  NT  <  '  D  <  I  >  -X 1  >  *S3  > 

920  Z2<I)=INT<<P<I>-Y1>*S4> 

930  OUTPUT  19  i"Pfl\Zl(I),Z2(I) 

;  " PD" 

940  NEXT  I 

950  OUTPUT  lO   "PUO, 980, 1O0, 90O 
it 

980  OUTPUT  10   "RIJ" 

970  OUTPUT  10  , "SI.22,  38" 

980  DISP  "INPUT  THE  LEGEND  YOU 

1990    INPL.T    p'iPE     'e'     WHEN    DGME"        A  label  or  legend  can  be  entered  here. 
2008     IF    P$="0"     THEN    2030 
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4.  FIBER  TEST  MODULE  (FTST) 

10     !  _FT::,T:::? 

2g     ;        JIM    NAGEOTTE  J  AN     13    87     The  Fiber  Test  is  known  as  "FTST". 

4  m    C L  E  A  R 

50    DISP    6    DISP    «YQU*MUST    FUN    Cfl  £*"*"  I2J!^IS?  b?- rU"  bCf°rC  lhC 
LB    BEFORE    PLOT    OR    SAM-DAT !  "     PLOT  or  SAV-DAT  routines. 

60    DISP    (2    DISP    "    USE    THE    FUNCTI 
ON    KEYS .  " 

70    SHORT    P (450 ) , L  ( 450 > , T 

SO  SHORT  D  1  '.  450  >  .  D2  (  450  >  ,  D3  (.  450 
>  •  D  <  4  5  O  ■' 

90  SHORT  Z 1  <  450  >  ,  Z2 '.  450  > 
100  GOTO  120 
1  10  CLEAR 

1  2  0    0  H    K  E  Y  #    4  ■  "  C  R  R  F  "    G  0  S  U  B    2  7  @  0 
1  3  O    0  N    K  E  V  #    1  ,  "  I  H  P  U  T  S  "    G  0  S  U  B    2  3  8 
140    0  H    K E  Y  #    2 ,  "  R C Q  U IP"    G  0 S  U B    6 O  O 
150    OH    KE^i  #    6."RD_DAT"    GOSUB     I  1  7  The  function  keys  are  defined. 

0 

1  6  O  6  N  K  E  Y  #  7  .  "  F  L  0  T  "  G  0  S  U  B  1  3  2  0 
170  ON  KEY#  5, "SflV_DPT"  GOSUB  29 

30 
1 8 O  ON  K E  Y ft  8 ."Cfl L B  "  G 0 S U B  133 8 
135  ON  KEY*  3, "AREA"  GOSUB  2500 

190    KEY    LABEL  The  program  loops  here  until  a  function  is 

200    GOTO    120  selected. 

210    CLEAR    @    DISP    "END    FTST87"    I? 

END 
220     ' 

230     !       Xt    INPUT':    ** 
240    CLEAR    e     !        DISP  "ENTER    DATA    R  The  "INPUTS"  routine  is  required  first. 

ATE<PTS  MIN>" 
250  !   INPUT  R 

260  R  =  6"00  !  PRESET  DATA  RATE 
270  DISP  "ENTER  ELAPSED  TINE  (30 

OP  4 5  >  " 
230  INPUT  EQ 
290  CLEAR  S  DISP  "ENTER  SPEED  <m 

m  rn  i  n  )  " 
300  INPUT  S 
310  CLEAR  G  DISP  "ENTER  GAGE  LE11  L  .    .  . 

r  T  H     (  mm  >  "  Thc  gage  lenSln  1S  t0  ^  uscd  in  computing 


320    INPUT    gi 

330  E1=E0*1000 

332  CLEAR 

350  DISP 

360  DISP  "ENTER  'A'  FROM  LDCALB 

LISTING" 
370  INPUT  C3 

330  DISP  "ENTER  ' B'  FROM  LDCALB" 
390  INPUT  B 
400  CLEAR  7^9 
410  OUTPUT  709  j "AI 10" 
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strain. 


4  20 

ENTEP  709  .;  '..' 

430 

DISP  "EXCITATION  VOLTAGE  IS" 

1  1 

4  40 

A=ABS<C3-  -'■.■> 

45U 

A=  I  NT  <  h 1 1  000  ...  -  1  000 

460. 

DISP   "h=":A 

4  70 

DISP  "B="-B 

475 

D  I  S  P  "  S  P  E  E  D  =  "  ;  S  :  "  M  M  /  M  I  N  " 

4  SO 

COPY 

496 

CLEAR  @  DISP  "ENTER  COMPLIAN 

CE  COEFPECIENTS" 

500 

DISP  "E=  " 

510 

INPUT  E 

5  2  O 

DISP  "F=  " 

530 

INPUT  F 

540 

DISP  "G=  " 

550 

INPUT  G 

560 

PRINT  "E=",E  i?  PRINT  "F="  ;F 

570 

PRINT  "G=";G 

zr  t-t  i-j 
_i  r,  M 

CLEAR 

530 

RETURN 

600 

i 

6  1  O 

!   **  AC  QUIP  %■% 

620 

IF  EO>0  THEN  668 

t>30 

BEEP  e  CLEAR 

640 

DISP  "'INPUTS'  MUST  BE  PUN  F 

IRST! " 

650 

RETURN 

6  6  0 

CLEAR 

690 

BEEP 

7O0 

FOP  1=1  TO  7  @  DISP  @  NEXT  I 

716 

DISP  "CHECK  IEEE  BUTTON.  MOU 

NT  SAMPLE  AND  BURN" 

720 

DISP  "HIT  CONT  WHEN  READY" 

722 

PAUSE 

723 

CLEAR 

724 

DISP  "FINDING  ZERO  LOAD" 

725 

OUTPUT  704  ■>  "K13,  "  ;  S 

730 

OUTPUT  709  .:  "ARVR1VT3VO5VA0A 

E2" 

OUTPUT  709  ;  "A  19" 

731 

732 

ENTER  709  ;  L 

733 

P=A*L+B 

734 

IF  P<0  THEN  GOTO  755 

735 

OUTPUT  704  :•   "K2" 

736 

WAIT  50O 

737 

OUTPUT  704  ;"K0" 

7  33 

GOTO  731 

750 

CLEAR  6  BEEP  S  DISP  "MOVING 

UP" 

755 

OUTPUT  704  : "K?" 

756 

OUTPUT  769  : "AIOH 

75? 

ENTER  709  ;  L 

759 

P=fl*L+B 

The  excitation  voltage  is  read  and  the 
calibration  adjusted  to  compensate  for  any 
voltage  drift. 


The  compliance  coefficients  are  entered  so 
that  a  matching  compliance  curve  may  be 
generated  and  then  subtracted  from  the  total 
deformation,  yielding  the  true  displacement 
of  the  sample. 


The  routine  ends  and  returns  to  the  menu. 

The  "ACQUIR"  routine  is  the  actual  data 
acquisition  routine. 

The  program  checks  to  see  that  the  inputs 
have  been  entered  and  provides  prompts  if  no 
input  data  can  be  found. 


The  OUTPUT  "704"  codes  provide  for 
automatic  control  of  the  Instron  tester. 
Initially  the  cross  head  is  adjusted  to  provide 
for  slight  tension  in  the  fiber,  after  which 
the  timer  will  start  and  data  collection  will 
begin. 
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760  IF  P<.06  THEN  GOTO  75b 

762  OUTPUT    704     , "K21"  The  voltage  from  the  load  cell  is  converted 

7  65  BEEP    >?    D  I  S  P     "TIMER    OH"  to  physical  units  of  grams.  Here  the  tension 
76ft  ON    TIMER*    1  ,  E 1    GOTO    830  level  is  adjusted. 

770  1=1 

7R0  OUTPUT  709   "RIO" 

790  ENTER  709  ;     L(I> 

800  1=1+1 

810  WHIT  23 

H2fi  GOTO  780 

830  OFF  TIMER*  1 

R35  nUTPUT  704  . "KO" 

840  1=1-1 _ 

8  50  WRIT  t-  O 
860  CLEAR 

870  FOR  T=l  TO  7  @  DISP  i?  NEXT  1 

880  BEEP  (3  CLEfiR  <*  DISP  "CHECK  E 

XTENSION     HIT  CONT" 

890  PR USE 

900  OUTPUT  704  ..  "Kl" 

910  N=I 

920  DISP  N,"ORTR  POINTS" 

930  DISP  "COMPUTING  " 

940  FOR  1=1  TO  N 

950  pc;i>  =  fi#L<ri>+B 

960  NEXT  I 

9  70  P2  =  P <  1  >  8  P 1  =P  <  1  ) 
980  FOR  1=2  TO  N 

990  IF  P< I)>P2  THEN  K=I 

1 0O0  IF  P < I > > P2  THEN  P2  =  P (  I  > 

1010  IF  PCIXP1  THEN  P1=P<I;- 

1030  NEXT     I 

1060  DISP    " P<max>  =  "  ;  P2  Aflcr  a  quick  analysis  of  some  key  data 

1070  D  4  =  S  *  (  K  '  R  >  *  <  I  • '  G  1  )  *  1 09  points  are  printed  out  to  provide  a  indication 

1  O  8  0  DISP    "  D ':  "  .;  D  4  .;  K  of  the  nature  of  the  data  gathered. 

1090  PRINT  "PCMflX>=  " ;P2 

1095  PRINT  "D":  =  "  ;D4,k 

1150  RETURN 

1160  ! 

1  170  "    **   RERD  DATA   :** 

1180  CLEAR 

1190  DISP  "ENTER  FILE  NAME" 

1195  DISP    "YOU    M  U  ST    HA  V  E    INPUT  S    The  sub  routine  returns  to  the  main  loop. 

BEFORE  YOU  CAN  PLOT  READ  DP 

T  A  "  The  "RD-DATA"  routine  allows  you  to  read 

1200  I  N  P  U  T    N  $  in  a  data  file  that  had  been  stored  earlier  and 

1210  ASSIGN*    1    TO    Nf  then  plot  results. 

1220  READ*    1  N.EO,S,Gi,k 

1230  FOP     1  =  1     TO    N  It  reads:  Number  of  data  points,        N 

1 249  READ*    1     ;     P  <  I  ;■ ,  D<  I  >  Elapsed  time  of  test,        EO 

1250  NEXT    I  Speed  of  test,  S 
IdbQ  ASSIGN*    I    tq    *                                                  Gage  length  tested,          Gl 
1270  r^=P<K.>    6    :;^  =  D(K)                                             Point  of  Max  load,           K 
1280  DISP    -    DATA    HAS    BEEN    READ"                   Load  and  Strain,             P(I),  D(I) 
1  c.  J  0  Bhbr 

13O0  RETURN 
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1 3  i  8  ! 

1320  '  1 1    PLOTTER    i'  *  Program  returns  to  the  menu. 

1330  CLEAR 

1340  DISP  "ENTER  Y  MAX  FOR  Y  AX  I 

S  "  The  "Plotter"  routine  is  the  same  as  used  in 

1350  INPUT    Y2  "MCTST". 

1360  DISP    "ENTER    V    MIN    FOR    Y    AX  I 

1370  INPUT  Yl 

1386  DISP  "ENTER  X  MAX  FOP  X  AX  I 

1330  INPUT  X2 

1400  DISP  "ENTER  X  MIN  FOR  X  AX  I 

14  10  INPUT  XI 

1420  S  3  = 1 0  O 0  < X 2 - X 1 > 

14  3  0  S  4  =  1  0  0  0  ' "'  Y  2  -  Y  1  ) 

1448  PRINTER  IS  10 

1450  CONTROL  18  5  •  4S 

14b0  OUTPU T  10  ; " I N , SP 1 ; I P2400 , 1 

b08  8868 • 6988 ; " 
1470  OUTPUT  10   "SCO, 1OO0, O, 10O0 

14S0  OUTPUT  18   "RNO, OPDO  ■  lOOO,  1 

800  •  1  O0O  ,  1  000  ■  O  •  8  ,  OPU  " 
1490  OUTPUT  10   "SIO  2, O  3, TL1 . 5 

,  O  " 

1  5  O  O  F  0  R  I  =  8  T  0  5 

1510  Y=INT<Y1  +  <Y2-Y1  >'5*I) 

1520  Y4=INT<<Y-'.  1  'tS4> 

1530  OUTPUT  18  , "PA  0, " , Y4 , " YT; " 

1540  OUTPUT  18  ; "CP-5, -O . 07; LB" ; 

Y  ;CHRf •  3  ■ 
1550  NEXT  I 
1560  FOR  1=1  TO  5 
1570  v  =  v 1  +  ( X2 - X !  )  / 5*1 
1 580  X4= I  NT  <  < X-: '  1 >  t  S3  > 
1590  OUTPUT  1 8   " P A " > X4 ,  " , 8 : XT ;  " 
16O0  OUTPUT  lO  , "CP-1 .3,-1; LB" ;X 

i CHR$<3> 
1610  NEXT  I 

1620  OUTPUT  10   "SI. 30.  42" 
1630  OUTPUT  10   "PA220.8;CP6,-2 

3  .:  LBSTRA  IN'  V.  >  "  ;  CHR$  <  3  > 
1648  OUTPUT  18  ,  "PAO, 460  =  010,  1; C 

P-2  6,2  60;  LB  FORCE  ';?rn>";C 

HR$<3> 
1650  OUTPUT  lO   "DI;PU" 
1660  FOR  1=1  TO  N 
1665  D-'I  >  =  S  *  I  ,R 

1680    Z2€i>*iNT<«P<i>-Vl>*l4>  Here:  D(I)  rePresents  *"  stra,n  whlle  lhe 

1690    OUTPUT    10     :  "PA"  .  Zl  <  I  >/22<  I  )  loadlsatP(I)- 
,  "PD" 
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1780  NEXT  I 

1716  OUTPUT  18  ;  "PU0,  900..  100,900 

ii 

1720  OUTPUT  10  .;  "PU" 

1730  OUTPUT  19  i "SI .22, .38" 

1748  DISP  @  OISP  "ENTER  THE  LEGE 

HO  YOU  WANT  " 
1758  DISP  "ENTER  A  'O'  WHEN  DONE 

II 

1760  INPUT  P* 

1770  IF  P*=,,8"  THEN  1800 

1  780  OUTPUT  10  .;  "  CP  .;  LB  "  .;  P$  ;  CHR*  ( 


Label  the  graph  for  identification. 


1790  GOTO  1""40 

1800  PRINTER  IS  2 

1810  RETURN 

1820  • 

18  30  '    *.**  CflLB  *** 

1840  IF  FO0  THEN  GOTO  1870 

1858  CLEAR  I?  DISP  "NEED  INPUTS'" 

I860  RETURN 

1870  CLEAR  @  BEEP  6  DISP  "DATA  f 

(INVERSION    TO    PHYSICAL    UNIT" 

i  &ft  q  FAR     1  =  1     T  fi    K  ^ne  rouune  returns  to  the  menu. 

1  890  0  I  <  I  ■  =  S *  I  sR  n„  a  r  nn  .    l       ,.       . 

19F10  D2  <L  I  >  =G  1  P  •'  I  '>  *-2  +  P  iP  <  I  >  +  E  CALB    is  the  calibration  routine  where  the 

1  9  1  0  D3  ■■  I  ■'  =0  1  (  I  ->  -02  <  I  >  ^'na'  outPul  IS  lTUC  displacement  and  strain 

1920  D<  I>=D3<  I  )<■- Gil  100     !     5S    STRfll    of  the  fiber  sample. 

N 

1930  NEXT     I  Dl   is  the  total  displacement.  D2  is  the 

1  940  FOP     I  =K+  1     T0    N  displacement  of  the  system  as  a  function  of 

195  0  0  1  '•  I  >  =  S  t  I     R  mc  l°ad,  PI,  calculated  using  the  coefficients 

I960  0  2  <  I  >  =  D  2  <  K  >  E,  F,  and  G. 

1970  03  <  I  >  =0  1  ■•  I  >  -02  CI)  D3  =  Dl  -D2,  the  true  displacement  of  the 

1988  0^1  >  =03  <  I  >  /  G 1  *  1 08     !     '-.    STR  A  I    fiber  sample. 

N  D  is  the  strain  in  percent . 

1990  NEXT    I 

2000  PRINT  "P(MflX>  =  " ;P<K)  ;K 

2010  PRINT  "05:  =".;0<K> 

2815  CLEAR  S  DISP  "CALB  COMPLETE 

li 

202O  RETURN 

2830  !  ***  STORE  ADJUSTED  DATA  * 

* 

204O  IF    O(30>O0    THEN    2870 

2850  BEEP    6    DISP     "    NO    DATA    PRESE    Return  to  the  menu. 

NT" 

2060  RETURN 

2878  CLEAR    8    BEEP    6    DISP    "INSERT   mis  sub  rouUne  is  used  to  store  die  adjusted 

2888  DISP    "ENTER    FILE    NAME."  data  calculated  in  the  CALB  sub  routine. 

2890  INPUT    N* 

2188  CLEAR  I?  DISP  "STORING  DATA" 

2118  M=2*<N+1 >*S+58 

2128  H=INT(M.'256>  +  1 

2138  CREATE  N$,H,256 

2148  ASSIGN*  1  TO  N$ 

2158  PRINT*  1  i     N,E8.S,G1,K 
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2160  FOP  1=1  TO  N 

2170  PRINT*  1  .;  PC  I  >,  DC  I) 

2130  NEXT  I 

2190  ASSIGN*  1  TO  * 

2200  CLEAR  @  DISP  "DATA  STORED" 

2210  BEEP 

2220  RETURN 

2500  ! 

25  10  !  **      AFEA     ** 

2520  CLEAR    I?    D  I  S  P     "COMPUTING    ARE  Return  to  menu 

A    UNDER    Px'D    CURVE  .  " 

2530  T  =  0  AREA  calculates  the  are  under  the  load 

2540  FOR     1  =  1     TO    K  -  1  versus  displacement  curve. 

2550  fll  =  CPCI>+PCI  +  l>>x'2 

2560  R2=fil*CDC  I  +  D-DCI)} 

2570  T=T+A2 

25S0  NEXT    I 

2590  PRINT     "AREA    UNDER    THE    P^D    C 

URVE="  i  T;  "  ■9m."'." 

2595  PRINT 

2688  IF  N*<>""  THEN  PRINT  ; " 

"  .  Nf 

2610  CLEAR  8  DISP  "AREA  COMPUTED 

"  @  RETURN 

2700  ! 

2710  '   tfDISP  CRT  GRAPH** 

27^0  PRINT     "  "  .=  Ht    »2      jhis  sub  routine  provides  for  a  graphical 

P  R  I  N  T  output  on  the  CRT  screen 

2730  CLEAR    @    DISP     "INPUT    PCMflX) 

FOP  Y  AXIS" 

2740  INPUT  P2 

2750  P1=0 

2760  GCLEAR 

2780  GRAPH 

2790  LDIP  0 

2800  SCALE  -43,203,-36,156 

2310  X AX  I  3  0,32 , 0 ,  1 60 

2820  VAX  IS  0,28,0.  140 

2830  FEN  1 

2340  PENUP 

2350  31  =  160 '-  C  E  O  •  6  O  >  /  S 

2860  S2=140/P2 

2870  FOP  T=l  TO  5 

2880  X=S*T'10 

2890  :  :  =  I  NT ( X* 1 00O  > / 1 00O 

2900  X2=X*S1 

2910  NOME  X2,-10 

2920  LABEL  VAL*O0 

2930  NEXT  T 

2940  FOR  T=0  TO  5 

2950  P  =  P2.-'5*T 

2960  P3=P*S2 

2970  MOVE  -23- P3 

2930  LABEL  VflLfCP) 

2990  NEXT  T 

3O00  MOVE  40,-23 

3810  LABEL  "DISP<MM>" 
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3820  DEG 

3S30  MO'-'E  -29/25 

3040  LDIP  90 

3050  LABEL  "FORCE  <*»>" 

3060  LDIR  0 

3U70  FOR  1=1  TO  N  ' 

3  0  8  0  D<I>=S*<I / R  > 

3090  ZKI>=D<I)*S1 

3100  Z2< I >=P< I )*S2 

3110  PLOT  Z1CI > , Z2CI) 

3120  NEXT  I 

3130  COPY 

3140  CLEAR 

3150  RETURN 

3200  !     **DIAMETER** 

2xi.:  ' .  DIAMETER    PROGRAM  This  sub  routine  is  used  when  dcicrmining 


fiber  diameter,  it  provides  for  the  input  of 
the  fiber  designation  and  micrometer 
readings. 

The  fiber  diameter  is  determined  based  on  an 


corrected  in  accordance  with  the  SEM  cal 
curve. 


3220  CLEAR 

3230  '  W=WflVE  LENGTH,  N*=FIBER  D 

ESI  GNAT  I  ON.  C1=DELTAC  AND  S 

3  =  SEM  EQUIVALENT  VALUE 
3240  OISP  "THIS  PROGRAM  FINDS  TH 

E  FINAL  DIAMETER" 

3250  W=    0000006328  average ^of  the  micrometer^ readings  and 

3251  N=l 

3252  C=115 

3253  L=698  1 

3254  DISP  "ENTER  FIBER  DESIGN" 
3260  INPUT  Nt 

3270  PRINT  "FIBER  NO    "iN$ 
3280  FOP  1=1  TO  H 

3290  DISP  "ENTER  MICRO  READING" 
3300  INPUT  CI 
33  1 0  Q=W*2*L-'  <  C  1  +C  ) 
3320  Q=IHTCQ*1  E 1 3 > / 16000096 
3330  PRINT  "HPSLT  '..'ALUE  =  " .;  Q 
3340  !  CONVERT  TO  SEM  VALUE 
3350  Q 1 = . 89XQ+ . 36 1 4 
336id    01  =  INT  <  01  *10000000>/ 1000000 

0 

3370  PRINT  "FIBER  DIAMETER=  ".;Q1 
; "  MICRO 

NS" 
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5.  LOAD  CELL  TEST  MODULE  (LDTST) 


1  0  i  M  LDTST  t.  i 

26  '         TAKE'?    £5    READINGS    FROM     LO  "LDTST"  is  a  short  utility  to  lest  ihc  load 
RDC E L L     RN  D    P  R I  N  T  S    fl '...' E  R  A  G  E  cell  for  proper  operation  and  reproducibility. 

30  L  L  E  A P 

40  DISP  "ENTER  h  " 

SO  INPUT  Rl 

SO  DISP  "ENTER  B" 

70  INPUT  Bl 

8  0  C  L  E  fl  R 

?  0  C  L  E  fl R  70? 

100  OUTPUT  70?  ;"AI10" 

1  10  ENTER  70?  ;  V 

I  20  R=fll'U 

125  B=B1 

130  P1=0 

1  4  ti  F  Li  R     I  =  1  ^  T  0    4'  5  Load  cell  output  voltage  is  recorded  25  times 

150  OUT  P  U  T    70?     ;  •'  R  1 0  "  then  converted  to  grams. 


The  average  load  is  printed. 


The  program  pauses  and  a  load  level  may  be 
changed. 


ISO  ENTER  70?  ,  L 

170  P=L*A+B 

1  75  P=I  NT  <  P  *  1  000  >  '  1  UOO 
ISO  DISP  P. 
1?0  F1=P1+P 

200  ne::t  i 

210  PRINT  "'vKout  >  =  "  ;  'J 

???  o  E  ^  r  T     "  L  ° R D    =  "  '  F  1 '"' 2 5  '■  "  ''  '9 "'  '  "     ^ss  continue  after  installing  new  load. 

4  3  y  B  E  E  P  6 

240  PRINT 

250  DISP  "SET  LORD.  HIT  CONT"    D    ,,.  ,  , 

.. , ..  _  Press    pause    or    reset    to  leave  program 

2  70  PRINT  l0°P- 
2  S  0  •  j  0  T  0    S  0 
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6.  LOAD  STORAGE  MODULE  (LOADSTOR) 


11 


9 
18 

11 
12 
15 

10 
10 
10 


T 
M 

! 

01 

CL- 
OU 

ON 
ON 
ON 
ON 
KE 
GO 
EN 

0  ! 

0  C 

0   0 


*  * 

LO 
RP 
CT 
J  I 
M 

S  U 
EA 
K 
K 
K 
K 
K 
Y 

TO 
D 

LE 
IS 


LOST  OP  :U 
PC'S  AND  STORES  DflTR  FILE 
ES  CREATED  WITH  FTST  AND 
ST 

M  NflGEOTTE   OCT  27  S6 
P 1  (.  4  5  0  )  .•  0  <  4  5  0  > 
B  5  0  0  0 
R 
FY*  1  .  "  P  D  -  C  "  G  0  S  U  B  1  0  0  0 


EY*  2 
EY#  3 
EY#  4 
EY#  8 
LABEL 

60 


"SA'v'-C"  GOSUB  2008 
"FO-F"  GOSUB  3000 
"SA'-'-F"  GOSUB  4001-3 
"END"  GOTO  150 


**  RD-C  ** 
AR  S  BEEP 
P  "  READ  COMPLIANCE  FILE 


1030  OISP 

1040  OISP  "INSERT  OATH  TAPE" 

1050  OISP  "ENTER  FILE  NAME" 

1060  INPUT  N* 

107O  CLEAR  &    DISP  "READING  DATA" 

1080  ASSIGN*  1  TO  N* 

1090  READ*  1  :  N,E0/S 

UO0  FOR  1  =  1  TO  N 

UIO  READ*  1  ;     PKI) 

1120  NEXT  I 

1130  ASSIGN*  1  TO  :*: 

1140  CLEAR 

1 150  DISP  "DATA  READ" 

1160  BEEP 

1170  RETURN 

2000  I  *X    SAVE  0  t* 

2010  CLEAR  @  DISP  "SAVE  COMPLIAN 

CE  FILE" 
2020  DISP 

2030  DISP  "INSERT  DATA  TAPE" 
2040  BEEP 

2050  DISP  "ENTER  FILE  NAME" 
2060  INPUT  Nf 
2070  CLEAR 

20S0  DISP  "SAVING  DATA" 
209O  M=<N+1 >*8+50 
2100  H=INT<M.-  256  >  +  ! 
2U0  CREATE  N$,H,25b 
2120  ASSIGN*  1  TO  N* 
2130  PRINT*  1  .;  N,EO..  S 
2140  FOR  1=1  TO  H 
2150  PRINT*  1  ;  PKI) 
2160  NEXT  I 
2170  ASSIGN*  1  TO  * 
2130  CLEAR 


"LOADSTOR"  is  a  utility  program  to  allow 
for  quick  storage  or  access  to  data  files.  This 
is  useful  for  the  transfer  of  files. 


'RD-C"  and  "SAV-C"  read  and  save  the 
'MCTST"  files. 


'RD-F"  and  "SAV-F"  read  and  save  the 
•FTST"  files. 


After  loading  a  data  file,  you  may  press 
"END"  and  then  examine  or  list  the  data  in 
the  command  mode.  Refer  to  the  print  out 
for  variable  names.  When  action  complete, 
re-enter  the  program  without  losing  the  data 
by  entering  "cont  50". 
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219@  DISP  "DATA  SAVED" 

220@  BEEP 

22  10  RETURN 

3600  i     **  RD-F  ** 

7010  CLEAR  i?  BEEP 

3020  DISP  "  READ  FIBER  TEST  FILE 

ii 

3  0  3  0  D  I  S  P 

3040  DISP  "INSERT  DATA  TAPE" 

3O50  DISP  "ENTEF  FILE  NAME" 

3060  INPUT  N* 

3070  CLEAR  S  DISP  "READING  DATA" 

30S0  ASSIGN*  1  TO  H$ 

309O  READ*  1  ;  N,E0,S,G1.K 

3100  FOR  1=1  TO  N 

3110  READ*  1  :  PKI),D(I) 

3120  NEXT  I 

3130  ASSIGN*  1  TQ  * 

3140  CLEAR 

3159  DISP  "DATA  READ" 

3168  BEEP 

3170  RETURN 

40OO  !  **  SAVE  F  ** 

401O  CLEAR  i?  DISP  "SAVE  FIBER  TE 

ST  FILE" 
402O  DISP 

403O  DISP  "INSERT  DATA  TAPE" 
4040  BEEP 

4O50  DISP  "ENTER  FILE  NAME" 
4O60  INPUT  N* 

4  070  C  L  E  A  R 

40SO  DISP  "SAVING  DATA" 

4090  M«2*<N+l)*S+50 

4  1O0  H=INT<M/256>+l 

4110  CREATE  N$.  H-256 

4120  ASSIGN*  1  TO  H$ 

4130  PRINT*  1  ;  N.E0,S.G1,K 

4140  FOR  1=1  TO  N 

4  150  PRINT*  1  ;  P1(I>#D<I> 

4160  NEXT  I 

4170  ASSIGN*  1  TO  * 

4180  CLEAR 

4  190  DISP  "DATA  SAVED" 

4200  BEEP 

4210  RETURN 

50OO  !  **  INFO  ** 

5010  PRINT  "APTER  'PAUSE'  OR  'EN 

D'  YOU  MAY  CONTINUE  WITHOUT 
DATA  LOSS  BY  TYPING = "       • 
5020  PRINT  " ' CONT  50 ' " 
5030  PRINT 
5040  PRINT  "C  FILES  CONTAIN  N , EO 

,S  AND  PKD" 
505O  PRINT 
5060  PRINT  "F  FILES  CONTAIN  N,EO 

, S,G1,K  AND  Pl<I)iO<D" 
5070  PRINT  @  PRINT 
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APPENDIX  C.  DETERMINING  BUNDLE  PRELOAD 

Before  composite  preload  level  can  be  determined  a  mathematical 
model,  based  on  individual  fiber  data,  must  be  developed  to  study  the  effects 
of  gage  length  on  the  strength  distribution.  Previous  studies  have  provided 
conclusive  results  that  composite  strength  is  a  function  of  gage  length.  In 
order  to  predict  the  behavior  of  a  254  millimeter  strand,  the  results  of  the 
fiber  tests,  (gage  length  =  50  millimeter)  must  be  examined  and  equated  to  a 
sample  that  is  254  millimeters  long. 

Results  from  a  previous  study  [  Ref.  7],  which  examined  the  affects  of 
proof  testing  Kevlar-49  at  elevated  temperatures  will  be  used  to  verify  the 
mathematical  model  and  to  establish  an  estimated  preload  level. 

1 .  MATHEMATICAL  CORRECTION  FOR  GAGE  LENGTH 

Fl  =  1 -exp{-(a/Pi)ail         for  gage  length  =  li  (4) 

Rl  =  exp{-o/Pi)«i  (5) 

12  /  li  =  n  for  gage  length  =  I2 

RT  =nRiR2-  =(Ri)n 
=  [exp{-a/p!)ai}}n 

=  exp  {-n(a/Pi)ai}        =   exp  { -  n(a /(32  )  °2  } 
=  exp  {  -n  (  aai  /  Piai )  }  =   exp  {  -n(aa2  /  (32a2)} 

acq  -  ca2    . .    ai  =  a2 

n/(3ia!  =  l/p2«2 

(Pl/p2)ai  =n  =  li/l2 

(Pl/fe)  =  l2(1/al>  /  liWai) 

(p2/pl)  =  (li/l2 )(!/«!) 


Results: 


P2    =   Pi    (  *1  /  h)  (1/al)  (6) 
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Mathematically,  equation  (6)  provides  for  the  beta  adjustment  for  gage 
length  . 


*&" 


2.      KEVLAR-49  MODEL 

Based  on  original  test  data,  Table  VIL,  the  beta  for  strain  pe  can  be 
corrected  for  various  gage  lengths. 

TABLE  VH  .    ORIGINAL  TEST  DATA 

Material  G.L.(mm)        a  P(gm)  (gm/£) 

Kevlar-49         50  6.7  38.9  1.3x103 

Kevlar(25oC)      254  25.9  9.74x103       4.0x105 

Kevlar  (70oC)       254  9.50x103        3.0x105 

Using  the  test  data  for  a  50  mm.  fiber  and  substituting  pe  for  pi  in 
equation  (6),  where; 

Pe  =   Jl 
E 

Pe2  =    Peldl/l2)(1/al) 
Correction  to  gage  lengths  can  be  made,  see  Table  VH[. 

TABLE  VIE.  CORRECTED  Pe  OF  KEVLAR  -  49  FOR 
VARIOUS  GAGE  LENGTHS   5. 

Delta  mm  5         .1         .5         1  10      50        254 

(ineffective  length) 

Beta  of  Strain  pe       .076    .059    .054    .038    .030     .024 

Using  the  reults  from  a  proof  test  of  a  Kevlar-49  strand  (gage  length  = 
254  mm)  at  70°C  the  number  of  individual  fiber  breaks  can  be  determined. 
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pe  =  Pl/£ 

Pe  =  9.50x103  gms  /  3.0x105  gms/e 

pe=  3.1x10-2 

Note  :  This  result  is  somewere  between  5=  10  mm.  and  50  mm.  (Table 
VIII.)  and  correcting  for  gage  length  ,  equation  (6),  S  =  50  mm. . 

The  percent  of  fiber  failures  F  as  a  function  of  a  possible  ineffective  gage 
length  is  determined  by  equation  (4). 

F=  l-exp{-(  (WPe70c    )aU 


TABLE  IX.  PERCENT  FIBER  BREAKS  AS  A  FUNCTION  OF 
OF  INEFFECTIVE  GAGE  LENGTH  70<>C 


5 

Pe 

F 

#of  fiber  breaks 
(267  per  bundle) 

.1 

.076 

.0001 

.03 

.5 

.059 

.005 

1 

1.0 

.054 

.01 

3 

10. 

.038 

.1 

27 

50 

.030 

.39 

104 

100 

.027 

.63 

168 

254 

.024 

.99 

267 

Based  on  the  fact  that  the  strand  was  heated  to  70°C  it  would  not  be 
unreasonable  to  expect  that  the  ineffective  length  5  approaches  50  mm. 
Therefore  from  Table  IX,  39  percent  of  the  fibers  broke. 
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3.      GRAPHITE  TESTING 

The  model  was  then  examined  for  graphite  fibers  based  on  known  test 
data  from  individual  fiber  testing  conducted  during  this  study.  For  gage 
length  =  50  mm; 

pL  =  15.5  gmf      a  =  4.0        pe  =  1.72  x  10-2mm/mm    E=  900  gm/e 

Using  equation  (6),  (3l  =  13.0gmf,  for  gage  length  =100  mm.  Now 
applying  this  information  to  a  dry  bundle  consisting  of  3000  fibers  254  mm. 
long  the  number  of  fiber  failures  can  be  predicted  as  a  function  of 
displacement  using  equation  (4),  (Table  X.) 

F=l  -exp{-(8/pe)ae} 
P£50mm  =  -0145  mm/mm 


TABLE  X.  PERCENT  FIBER  BREAKS  AS  A  FUNCTION  OF 
STRAIN  FOR  GAGE  LENGTH  OF  254  MM. 

£  AL  (mm)                    F  No.  of  Fiber  Failures 

.001  .254  .00006  <1 

.002  .508  .0009  3 

.004  1.016  .015  44 

.006  1.524  .07  214 

.008  2.032  .21  626 

.010  2.54  .46  1306 

.012  3.04  .69  2083 

.014  3.56  .89  2666 

.016  4.06  .98  2929 

.018  4.57  .99  2992 

.020  5.08  .999  2999 


For  example,  based  on  what  is  know  of  the  successful  Kevlar  test  at  70 
°C  (5=50  mm)  it  is  desired  to  break  40  percent  of  the  fibers  during 
preloading.  So  now,  the  graphite  bundle  would  be  loaded  until  delta  length 
equals  2.5  mm.  (Table  X). 
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These  results  can  be  graphical  confirmed  by  reviewing  a  representative 
load-displacement  plot  (Figure  18)  of  a  dry  bundle  (gage  length  =  254  mm.). 


ID 

(& 

o 


Apporx.    40% 
fibers  broken 


1,0  2.0  3.0 

Displacement  mm 


4,0 


Figure  18       Load-Displacement  Plot 
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APPENDIX  D.  TABULATED  TEST  RESULTS 


TABLE  XI. 

SINGLE  GRAPHITE  FIBER 

DIAMETER  RESULTS  (MICRONS) 

008  SERIES 

019  SERIES 

6.707 

7.282 

6.208 

7.291 

6.784 

7.283 

6.435 

7.300 

6.797 

7.285 

6.515 

7.310 

6.800 

7.285 

6.593 

7.314 

6.815 

7.290 

6.730 

7.316 

6.852 

7.294 

6.798 

7.333 

6.866 

7.306 

6.804 

7.335 

6.898 

7.307 

6.815 

7.338 

6.918 

7.318 

6.825 

7.345 

6.940 

7.321 

6.889 

7.350 

6.964 

7.322 

6.915 

7.350 

6.965 

7.324 

6.945 

7.351 

6.973 

7.330 

6.962 

7.364 

6.983 

7.335 

6.996 

7.372 

6.997 

7.336 

6.998 

7.372 

6.997 

7.337 

7.015 

7.380 

7.020 

7.341 

7.030 

7.389 

7.039 

7.358 

7.042 

7.395 

7.045 

7.359 

7.054 

7.426 

7.060 

7.370 

7.062 

7.430 

7.066 

7.376 

7.062 

7.440 

7.084 

7.377 

7.070 

7.447 

7.104 

7.377 

7.101 

7.453 

7.115 

7.378 

7.137 

7.454 

7.140 

7.388 

7.145 

7.458 

7.141 

7.397 

7.147 

7.480 

7.142 

7.403 

7.178 

7.482 

7.179 

7.420 

7.180 

7.507 

7.184 

7.425 

7.197 

7.516 

7.209 

7.445 

7.199 

7.551 

7.227 

7.456 

7.204 

7.569 
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7.236 
7.241 
7.243 
7.255 
7.258 
7.270 
7.275 


7.473 

7.206 

7.585 

7.573 

7.208 

7.586 

7.576 

7.212 

7.589 

7.650 

7.212 

7.600 

7.682 

7.216 

7.623 

7.760 

7.219 

7.670 

7.836 

7.223 

7.738 

7.227 

7.792 

7.243 

7.825 

7.289 

7.871 
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TABLE  XII.    SINGLE  GRAPHITE  FIBER 
FAILURE  LOAD  RESULTS 

019-N  019-1  008-N  008-1 


8.096 

4.303 

10.592 

8.022 

9.284 

5.015 

10.654 

9.255 

9.734 

6.819 

10.776 

9.653 

10.219 

7.368 

11.451 

10.707 

10.342 

8.111 

11.483 

11.018 

10.772 

9.899 

12.043 

11.157 

11.686 

10.324 

12.123 

11.243 

12.505 

10.636 

12.141 

11.622 

12.690 

11.751 

12.252 

11.781 

12.768 

11.770 

12.398 

11.804 

13.180 

12.440 

12.439 

12.499 

13.196 

12.448 

12.664 

12.623 

13.457 

12.667 

12.858 

13.071 

13.762 

12.904 

13.604 

13.157 

13.872 

13.703 

14.440 

13.297 

13.980 

13.705 

14.614 

13.630 

14.271 

14.265 

15.107 

14.140 

14.528 

14.502 

15.209 

14.812 

14.678 

14.528 

15.241 

14.914 

14.767 

14.811 

15.372 

14.936 

15.033 

15.565 

15.567 

16.116 

15.521 

16.957 

15.851 

16.369 

15.522 

17.073 

16.270 

16.633 

15.615 

17.726 

16.272 

16.941 

15.956 

17.804 

16.664 

17.560 

16.145 

18.259 

17.352 

17.607 

16.242 

18.801 

17.360 

17.923 

16.369 

18.956 

18.009 

18.992 

17.545 

19.392 

18.064 

20.603 

18.081 

19.556 

18.586 

20.911 

19.962 

19.610 

19.544 

21.572 

20.352 

19.915 

20.490 

25.442 

20.627 

N  =  NPSIT  Results  I  =  INSTRON  Results 
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TABLE 

Xm.  FAILURE  LOAD  RESULTS  OF 
SERIES  019  GRAPHITE 

TRIAL 

LOAD 

DIAMETER 

TRIAL 

LOAD 

DIAMETER 

(I) 

(gms) 

(microns) 

(N) 

(gms) 

(microns) 

1 

3.447 

7.219 

1 

7.284 

7.585 

2 

4.303 

7.440 

2 

8.096 

6.815 

3 

5.015 

7.426 

3 

9.284 

6.800 

4 

6.819 

7.430 

4 

9.734 

7.242 

5 

7.368 

6.992 

5 

10.219 

7.350 

6 

8.111 

6.998 

6 

10.342 

6.515 

7 

9.899 

7.364 

7 

10.772 

7.344 

8 

10.324 

7.145 

8 

11.686 

6.593 

9 

10.636 

7.054 

9 

12.505 

7.227 

10 

11.751 

7.335 

10 

12.690 

7.482 

11 

11.770 

7.333 

11 

12.768 

7.516 

12 

12.440 

6.730 

12 

13.180 

7.589 

13 

12.448 

7.030 

13 

13.196 

7.458 

14 

12.667 

7.042 

14 

13.457 

7.453 

15 

12.904 

7.147 

15 

13.762 

7.507 

16 

13.703 

6.804 

16 

13.872 

7.600 

17 

13.705 

7.289 

17 

13.980 

7.447 

18 

14.265 

7.197 

18 

14.271 

7.015 

19 

14.502 

6.825 

19 

14.528 

7.454 

20 

14.528 

7.825 

20 

14.678 

7.380 

21 

14.811 

7.569 

21 

14.767 

7.206 

22 

15.565 

6.945 

22 

15.003 

7.551 

23 

16.957 

7.062 

23 

15.521 

7.137 

24 

17.073 

7.310 

24 

15.522 

6.889 

25 

17.726 

7.291 

25 

15.615 

7.350 

26 

17.804 

7.316 

26 

15.596 

7.623 

27 

18.259 

7.738 

27 

16.145 

7.101 

28 

18.801 

7.351 

28 

16.242 

7.208 

29 

18.956 

7.372 

29 

16.369 

7.872 

30 

19.392 

7.199 

30 

17.545 

7.212 

31 

19.556 

7.216 

31 

18.081 

7.338 

32 

19.610 

7.070 

32 

19.962 

7.389 

33 

19.915 

6.962 

33 

20.352 

7.586 

34 

20.627 

7.792 
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TABLE 

XIV.  FAILURE  LOAD  RESULTS  OF 
SERIES  008  GRAPHITE 

TRIAL 

LOAD 

DIAMETER 

TRIAL 

LOAD 

DIAMETER 

(D 

(gms) 

(microns) 

(N) 

(gms) 

(microns) 

1 

5.504 

7.283 

1 

10.234 

6.832 

2 

8.022 

7.388 

2 

10.592 

7.336 

3 

9.255 

6.852 

3 

10.654 

6.707 

4 

9.356 

6.815 

4 

10.776 

7.403 

5 

9.653 

6.965 

5 

11.451 

6.964 

6 

10.707 

7.241 

6 

11.483 

7.445 

7 

11.018 

6.918 

7 

12.043 

7.236 

8 

11.157 

7.066 

8 

12.123 

7.140 

9 

11.243 

7.285 

9 

12.141 

7.039 

10 

11.622 

7.084 

10 

12.252 

7.321 

11 

11.781 

7.060 

11 

12.398 

7.282 

12 

11.804 

7.330 

12 

12.439 

7.318 

13 

12.499 

7.377 

13 

12.664 

7.141 

14 

12.623 

7.576 

14 

12.858 

7.307 

15 

13.071 

7.358 

15 

13.604 

7.337 

16 

13.157 

7.306 

16 

14.440 

7.322 

17 

13.297 

7.285 

17 

14.614 

7.275 

18 

13.630 

7.258 

18 

15.107 

6.784 

19 

14.140 

7.224 

19 

15.209 

7.294 

20 

14.812 

7.184 

20 

15.241 

7.378 

21 

14.914 

7.045 

21 

15.372 

7.115 

22 

14.963 

7.142 

22 

15.567 

7.376 

23 

16.116 

6.973 

23 

15.851 

7.270 

24 

16.369 

7.020 

24 

16.270 

7.425 

25 

16.633 

6.944 

25 

16.272 

7.179 

26 

16.941 

7.456 

26 

16.664 

7.377 

27 

17.560 

7.397 

27 

17.352 

7.359 

28 

17.607 

7.473 

28 

17.360 

7.209 

29 

17.923 

7.335 

29 

17.360 

7.243 

30 

18.992 

7.420 

30 

18.009 

7.836 

31 

20.603 

6.997 

31 

18.064 

6.997 

32 

20.911 

7.682 

32 

18.586 

7.324 

33 

21.572 

7.642 

33 

19.544 

6.800 

34 

25.442 

7.760 

34 

20.490 

6.983 
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TABLE  XV.    COMPOSITE  STRAND  FAILURE 
LOAD  RESULTS 


PRELOADED 

(KG) 

43.17 

46.76 

43.38 

46.79 

43.41 

46.95 

43.81 

47.28 

44.40 

47.46 

44.75 

47.84 

44.94 

47.97 

45.69 

49.02 

45.72 

49.91 

45.74 

49.99 

46.04 

50.17 

46.42 

sfON-PRELOADED  * 

(kg; 

i 

40.80 

» 

47.29 

42.18 

47.54 

42.38 

47.69 

43.37 

48.05 

44.48 

48.31 

44.86 

48.44 

44.89 

49.08 

45.18 

50.35 

45.61 

50.60 

46.53 

50.80 

46.79 

51.91 

47.15 

52.22 

*  Preliminary  data  on  AS4-Graphite  Strand  Intrinsic  Strength, 

By:  E.M.  Wu,  N.Q.  Nguyen  and  G.W.  Nypiuk,  Lawrence  Livermore 
National  Labratory,    Livermore,  California,  April,  1972 
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